Statistical signatures of critical behavior in small systems 
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The cluster distributions of different systems are examined 
to search for signatures of a continuous phase transition. In 
a system known to possess such a phase transition, both sen- 
sitive and insensitive signatures are present; while in systems 
known not to possess such a phase transition, only insensitive 
signatures are present. It is shown that nuclear multifragmen- 
tation results in cluster distributions belonging to the former 
category, suggesting that the fragments are the result of a 
continuous phase transition. 

PACS number(s): 25.70 Pq, 64.60. Ak, 24.60.Ky, 05.70. Jk 



I. INTRODUCTION 

Beginning in the 1970's significant advances in the 
understanding of nuclear multifragmentation were made 
possible with the advent of high statistics inclusive ex- 
periments. Typically, only one intermediate mass frag- 
ment ( 3 < Zf < 30 ) was detected per event. From 
these inclusive studies came the first evidence that in- 
termediate mass fragments (IMFs) were associated with 
a simultaneous multi-body breakup of a system which 
had undergone expansion. A study of the fragment mass 
yield distribution obtained in an inclusive gas jet exper- 
iment conducted at Fermilab contained the first indica- 
tion that nuclear multifragmentation might be related 
to critical phenomena normally observed in macroscopic 
systems §. The Purdue Group was the first to make the 
suggestion that the observed power law in the fragment 
yield distribution might result from a system whose exci- 
tation energy was comparable to its total binding energy 
0. The exponent of the power law was 2<r<3, within 
the range expected for a system near its critical point. 



The presence of the power law and the value of the ex- 
ponent, coupled with the strong similarity of the nuclear 
and van der Waals potentials, led the Purdue group to 
suggest that multifragmentation of nuclei might be anal- 
ogous to a fluid undergoing a continuous phase transition 
from a liquid to a gas. Furthermore, the Fisher Droplet 
Model (FMD) 0], used to describe condensation in 
a fluid system near its critical point, after modification 
for nuclear physics effects, was capable of describing the 
isotopic yields of 50 fragments with one set of parameters 
0,0- The temperature of the system was determined 
to be about 5 MeV 0, a reasonable value considering 
that the average binding energy per nucleon in a nucleus 
is approximately 8 MeV. The success of this approach 
reinforced the notion that multifragmentation was both 
a thermal process and that it was related to critical phe- 
nomena. 

With the advent of exclusive experiments capable of 
detecting all of the charged reaction products, the pos- 
sibility of studying multifragmentation on an event-by- 
event basis became a reality. High statistics exclusive 
experiments in which the fragmenting system is charac- 
terized according to its nucleon number and excitation 
energy permit both the correlation of dynamical and sta- 
tistical information and the study of fluctuations in ex- 
perimental observables. Fluctuations are central to all 
critical phenomena, and indeed, such fluctuations are ap- 
parent in exclusive multifragmentation data. In this pa- 
per, the focus will be on the statistical signals of multi- 
fragmentation data observed in the EOS experiment 
]l0| . Comparisons will be made with two other systems, 
one of which exhibits critical behavior and one of which 
does not. 
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Much of the pioneering work in understanding the sta- 
tistical aspects of multifragmentation has been performed 
by Campi 0- @ and Mekjian @- ||. Both ef- 
forts have compared multifragmentation data to model 
systems in order to gain some insight into the nuclear 
breakup process. In this paper, many of the ideas sug- 
gested by these authors are followed and applied to both 
the EOS data and the model systems in order to demon- 
strate which of the many suggested signals are useful for 
the identification of critical behavior. A major goal of 
this paper is to present a comprehensive review of sev- 
eral methods proposed for detecting signals of critical 
phenomena in multifragmentation. 

It is tempting to compare the experimental data to dy- 
namical models that attempt to describe nuclear multi- 
fragmentation. However, the task of modeling multifrag- 
mentation from the initial collision phase of the reaction 
to freeze-out has proven to be a daunting task. Models 
that adequately describe the initial stage of the reaction 
p6|- HH do not satisfactorily describe the fragment for- 
mation stage, in either statistical or dynamical aspects. 
Likewise, the most successful models in describing the 
statistical properties of nuclear multifragmentation ]32|]- 
pTj] , assume thermodynamic equilibrium, yet fail to ad- 
equately match the dynamical features of the data. 

Molecular dynamical approaches, which have enjoyed 
considerable success in describing critical behavior in 
classical systems [t 



1 40 



have not been conclusive in 
describing nuclear multifragmentation and at times have 
yielded contradictory results Q , |4l]] . Later studies sug- 
gested flaws in the application of molecular dynamical 
models to nuclear multifragmentation, therefore calling 
into question the conclusions drawn from the earlier stud- 
ies (§. 

The most striking of the early theoretical efforts came 
from Campi's analysis of a few hundred compl etely re- 
constructed emulsion multifragmentation events |43[ and 
the comparison of these data to clusters generated from 
a percolation calculation (ll), [jl2| . In this series of pa- 
pers it was shown that the fragment distributions from 
multifragmentation bore a striking similarity to the clus- 
ter distributions from percolation lattices. This analysis 
provided strong evidence that multifragmentation was a 
statistical process which appeared to be related to critical 
phenomena. In this analysis another estimate of the ex- 
ponent t was made which agreed with the first measure- 
ments from the Purdue Group and several later analyses 
of various fragment distributions. 

In the early 1990's the ALADIN Group from GSI per- 
formed several multifragmentation experiments Q]- fl46| . 
Of particular importance was the "rise and fall" of multi- 
fragmentation. In one analysis the ALADIN group plot- 
ted the "rise and fall" curve of the production of IMFs 
versus an observable related to the excitation energy of 
the reaction for several multifragmenting systems. With 
the appropriate scaling the data collapsed to a single 



curve suggesting that the multifragmenting systems re- 
tained no memory of the reaction entrance channel. This 
is expected for an equilibrated system. 

The results of some statistical analyses of multifrag- 
mentation data could be interpreted to suggest that mul- 
tifragmentation is a sequential decay [E^ in contrast to 
the phase transition picture. The same sort of statistical 
analysis has also been applied to explicitly simultane- 
ous models Q and produced results that were similar 
to those of multifragmentation data. Thus those signals 
could be interpreted as evidence for either sequential or 
simultaneous multifragmentation p7[ . 

This last effort puts into focus the main question in 
this work: what type of analysis of the statistical aspects 
of a cluster distribution can provide the most insight into 
the nature of the mechanism which created the clusters? 
Specifically, can those systems which contain critical be- 
havior be distinguished from those which do not? It will 
be argued that this question has two answers. Analy- 
sis of the insensitive features of the cluster distribution 
cannot make the above mentioned distinction (i^] . How- 
ever, an analysis of the sensitive features of the cluster 
distribution will be shown to provide deeper insight into 
the cluster production mechanism. This type of analysis 
has been previously reported for clusters resulting from 
nuclear multifragmentation ||- j[o|. Note that the more 
generic term cluster will be used to refer to any compos- 
ite of constituents, whether these be molecules of a fluid, 
nuclear fragments or percolation clusters. 

The method employed to address the above question is 
as follows. The same analysis is performed on the clus- 
ter distributions produced by three different systems. In 
one case, clusters are generated by randomly partitioning 
an integer. Such one-dimensional partitioning does not 
posses critical behavior indicative of a continuous phase 
transition. In the second case, three-dimensional bond 
building percolation is used to produce clusters. Percola- 
tion is well-known mathematical construct that possesses 
a continuous phase transition, i.e. a critical point. Fi- 
nally, the cluster distributions resulting from the multi- 
fragmentation of gold nuclei arc analyzed. Although it is 
not known, a priori, whether the nuclear multifragmen- 
tation bears any relation to critical phenomena, it will 
be seen that the analysis presented in this work yields 
suggestive results. 

This paper is organized as follows. In section II a 
description of each system is presented. In section III 
the Fisher Droplet Model is reviewed. In section IV- 
A the insensitive signatures of the cluster distributions 
for all systems are examined. In section IV-B the sensi- 
tive signatures are examined. Sections V and VI present 
possible corrections to the analysis of the multifragmen- 
tation data. Finally, Section VII discusses the conclu- 
sions reached upon the completion of the analyses in sec- 
tions IV and V. Throughout this paper the term continu- 
ous phase transition will be used instead of second order 
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phase transition, the latter from the outdated Ehrenfest 
theory of phase transitions jf9| . 

II. DESCRIPTION OF SYSTEMS UNDER STUDY 
A. 1.0 A GeV Au + C multifragmentation 

Approximately 40,000 fully reconstructed events (76 < 
Zobserved < 82) were collected with the EOS experimen- 
tal apparatus discussed in ref . Q . In the collision of the 
projectile gold nucleus (197, 79) and the target carbon 
nucleus, so-called prompt nucleons are knocked out of the 
gold nucleus by quasi-elastic and inelastic collisions be- 
tween projectile and target nucleons ]5C[| . Immdeiately 
following the collision, he gold projectile remnant is in 
an excited state with fewer than 197 nucleons. The ex- 
cited remnant cools and expands, evolving to the neigh- 
borhood of the critical point in the temperature-density 
plane |5l| , where clusters condense from a high tempera- 
ture low density vapor of nucleons. The charge and mass 
of the projectile remnant, Zq and Aq, were determined 
for each event by subtracting the charge and mass of the 
prompt particles from the charge and mass of the gold 
nucleus JsT]]. Prompt particles have Zf = 0, 1 and 2 and 
are removed from the cluster distributions analyzed in 
this work. Only clusters created from the excited gold 
projectile are considered in the ensuing analysis. For 
events with the lowest charged particle multiplicities, to, 
the remnant had Z ~ 76, Aq ~ 194 and E* /Aq — 2 
MeV/nucleon, while for events with the highest multiplic- 
ities the remnant had Zq ~ 39, Aq ~ 92 and E* j Aq ~ 16 
MeV/nucleon 

Clusters of a given charge, Zf, were counted on an 
event by event basis to determine the cluster charge dis- 
tribution, Nzf ■ In this analysis, although the mass num- 
ber of the clusters is of interest, a cluster's charge will be 
used as an index. Mass numbers for clusters of charge 
one and two were measured in the EOS time projec- 
tion chamber. Clusters with Zf > 3 were assigned a 
mass number, Af, by multiplying the cluster charge by 
the mass to charge ratio of the excited gold projectile 
remnant; for low to events Aq/Zq ~ 2.55 and for high 
to events Aq/Zq ~ 2.36. This procedure provided an 
estimate of a cluster's mass number prior to any sec- 
ondary decay effects. It is assumed that on average 
Naj — Nzf ■ Finally, it is the normalized cluster distri- 
bution, nA f = -/Va,/j4o(to), that is used in the analysis 
presented in this paper. 

B. Percolation 

Bond building percolation calculations were performed 
on three dimensional simple cubic lattices of 216 sites. 
Cluster distributions for 100, 000 lattice realizations were 



generated in the standard fashion by forming bonds be- 
tween sites. Bonds were either active (on) or inactive 
(off) according to the following algorithm. 

The control parameter (e.g. temperature in thermody- 
namic systems) for percolation is the lattice probability, 
pi. A single value of pi was chosen for the entire lat- 
tice. All probabilities were between and 1. Next, a 
bond probability, pj, ( , was randomly chosen from a uni- 
form distribution on (0,1) for the i th bond. If p\ )i was 
less than pi, then the i th bond was active and two sites 
were joined into a cluster. This process was performed 
for each bond in the lattice. 

At low values of pi, few bonds were formed resulting 
in a high multiplicity, to, of small clusters, a distribution 
analogous to the gaseous phase of a fluid. At high values 
of pi , many bonds were formed resulting in a low mul- 
tiplicity of mostly large clusters, analogous to the liquid 
phase of a fluid. In an infinite lattice the phase transi- 
tion occurs at a unique value of the lattice probability, 
p c , when the probability of forming a percolating cluster 
changes from zero to unity. 

To examine the behavior of the average cluster distri- 
bution, the number of clusters of size A f per lattice site 
was calculated by histogramming the 100, 000 lattice re- 
alizations into 100 bins from to 1. The use of to as a 
control parameter and the ensuing effects on signatures 
of continuous phase transition were investigated by cal- 
culating the average number of clusters of size Af with 
the 100, 000 lattice realizations histogrammed in units of 

TO. 

C. Random partitions 

Random partitions were generated from 79 total sys- 
tem constituents, chosen to approximate the number of 
charges in the gold multifragmentation system. The al- 
gorithm is as follows. First a random choice of to was 
made from a uniform distribution on (1,79). Next the 
maximum size of a cluster, A^^, for an event with to 
was determined; this depended on the constraints of the 
system size, Aq = 79 and the choice of to. The size of the 
first cluster, A±, was then randomly chosen from a uni- 
form distribution on (1, -4*, oa .). There were then to — 1 
clusters to be generated from 79 — Ai constituents. The 
maximum size of a cluster for an to — 1 event from a 
79 — 1 constituent system was determined: A^ ax . The 
size of the second cluster, A2, was then randomly chosen 
form a uniform distribution on (l,A^ aa .). This process 
was repeated until all constituents belonged to a clus- 
ter. 100, 000 partitions were generated in this manner. 
This particular weighting results in a power law cluster 
distribution Q. 
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III. REVIEW OF THE FISHER DROPLET 
MODEL 

The focus of most studies of phase transitions is on 
standard thermodynamical variables such as a system's 
temperature, density, compressibility, etc. These quan- 
tities are difficult or impossible to measure directly in 
present nuclear multifragmentation experiments. Thus 
a theory which addresses quantities accessible to MF ex- 
periments is needed. To that end Fisher's gas-to-liquid 
phase transition model, based on Mayer's condensation 
theory, is followed §), g, |§. 

Fisher begins his model, called the Fisher Droplet 
Model (FDM) hereafter, by writing the free energy for 
the formation of clusters of size Af as: 

AG Af = -k b TA f \n(g{(x,T)) - 

k b T]n(f(A f ,T)) + k b Trhx{A f ) + ... (1) 

Where k b is the Boltzmann constant and the g-term is 
the bulk formation energy, or volume term and: 

g(p,,T) = exp[(fi - Hcoex)/k b T], (2) 

where [i is the chemical potential and p. CO ex is the chem- 
ical potential along the coexistence curve. 

The /-term is related to the surface free energy of clus- 
ter formation. It's a form given by Fisher is: 

f(A f ,T) = exp[a u>A'}eT c /k b T], (3) 

where a is a critical exponent and is related to the ratio 
of the dimensionality of the surface to the dimensionality 
of the volume, ao is a constant of proportionality relat- 
ing the average surface area of a droplet to its number 
of constituents and u> is the surface entropy density; e is 
a measure of the distance from the critical point. For 
usual thermodynamic systems e = (T c — T)/T Cl in the 
percolation treatment e = (pi —p c )/Pc and for multifrag- 
mentation e = (m c — m) / m c will be used. All formula- 
tions of e are such that e > (e < 0) corresponds to the 
liquid (gas) region. This form of the surface free energy is 
applicable on only one side of the critical point, the sin- 
gle phase side. A more general form suggested by efforts 
from percolation theory J53|- |36| that can be applied on 
both sides of the critical point and leads to a power law 
which describes the behavior of the order parameter is: 

f(z) = Aexp[-(z-B) 2 /C], (4) 

where the scaling variable, z, is 

z = A a f e. (5) 

The physical interpretation of the parameters A, B and 
C is an open question. 

Finally r is another critical exponent depending prin- 
cipally on the dimensionality of the system and has its 



origins in considerations of a three dimensional random 
walk of a surface closing in on itself. For three dimen- 
sions 2 < t < 3 J57|. In eq. (|l|), qg is a normalization 
constant which will be shown to depend solely on the 
value of t H). 

From the free energy of cluster formation the average 
cluster distribution normalized to the size of the system 
is: 

n Af (e) = exp(~AG Af /k b T) = q Aj T f(z)g(p,T) A i . (6) 

At the critical point, e = 0, both / and g are unity and 
the cluster distribution is given by a pure power law: 

n A f {z) = qoAj T . (7) 

If the first moment of the normalized cluster distribution 
is considered at the critical point then [fi8|| : 

Mi ( e = o) = Y, Af n *f(t) A f = ioj2 Af Al r = lo ( g ) 

when the sum runs over all clusters. From eq. (^) it 
is obvious that the value of the overall cluster distribu- 
tion normalization constant, qo, is dependent on r via a 
Riemann ^-function: 

The above is true only if the scaling assumptions in the 
FDM apply to all clusters. For finite size systems even at 
the critical point this is only approximately true. How- 
ever, it will be seen that eq. holds reasonably well 
at the critical point for systems with a continuous phase 
transition over some range in cluster size. 

In the FDM it is assumed that all clusters of size Af 
can be treated as an ideal gas, so that the total pressure 
of the entire cluster distribution can be determined by 
summing all of the partial pressures: 

P/(k b T) = Y JA n Af {e) 

= <loJ2 Af A } T f( z )9(P> T ) Af = M o(e) (10) 

Is is clear from eq. (|l(]) that the pressure of the system is 
related to the zeroth moment of the cluster distribution. 
The density is then: 

= Y. Af nA ^ A f = M ^- ( n ) 

The density is given by the first moment of the cluster 
distribution. 

It is now a simple matter to derive the power law which 
describes the divergence of the isothermal compressibil- 
ity, kt- By definition: 



4 



18V I dp 

KT = -v { dP )T = -p { dP )T - 

Noting that k b Tp = g{p,T){dP/dg(p,T)), eq. 
be rewritten as: 



(12) 
21) can 



-1 



kt = 



9(M,T)- 



dP 



g{p,T) 



d 2 p 



■dg(p,T) ^' ' dg(p,Tf 
which leads to: 

kt = (phT)- 1 + {p 2 k b T)- 1 V n Af (e)A 2 
= {pk b T)~ l + {p 2 k b T)- 1 M 2 {e). 



(13) 



(14) 



The sum in the second term illustrates the relation of the 
second moment of the cluster distribution, M 2 (e), to the 
isothermal compressibility. The sums in eq. (|To|), ( pd| ) 
and ( fi~4] ) run over all clusters in the gas and exclude the 
bulk liquid drop. In percolation and multifragmentation 
the largest cluster on the liquid side of the critical point 
will be considered as the liquid drop and will thus be 
excluded from the sum. On the gas side of the critical 
point, the sum runs over all clusters as there is no longer 
a liquid drop. 

In the thermodynamic limit, large Af dominate the 
sum so that it may be treated as an integral giving: 



k t = (phT)- 1 + (p 2 k b TY 



n Af (e)A 2 f dA f . (15) 



Working along the liquid-gas coexistence curve so that 
g(p,T) = 1 eq. ( |l5|) reduces to: 



kt 



(phT)- 1 + (p'hT)- 1 



A 2 f - T f(z)dA f . (16) 



A change of variables from Af to z shows that near the 
critical point: 



Kt 



[p 2 k b T)- x 



qo 

° JO 



dz f(z) 



(17) 



This is the so-called 7-power law which describes the di- 
vergence of the isothermal compressibility and the second 
moment of the cluster distribution near the critical point. 
The scaling relation between the exponents 7, a and r 
is: 



7 = 



(18) 



The absolute normalization constants of the Mi{t) 
power law depend on the scaling function, f(z), the ex- 
ponent a and the overall normalization of the cluster dis- 
tribution, go, which in turn depends on the exponent r: 



r+ = 



±00 



dz f(z) \z\ 



(19) 



The second moment is related to the isothermal com- 
pressibility by the temperature and density of the sys- 
tem. 

The derivation of the 7-power law demonstrates one 
way to arrive at the scaling relations between the critical 
exponents. In addition it illustrates the existence of only 
two independent exponents and shows the relation of the 
moments of the cluster distribution to familiar thermo- 
dynamic quantities. Fisher's framework here illustrated 
and tempered by percolation theory will be used in the 
analysis of the cluster distributions of the three systems 
discussed above. It will be seen that in the case of sys- 
tems which exhibit a continuous phase transition, the 
framework of Fisher is well followed, while for systems 
with no such phase transition, the framework fails, as it 
should. 



IV. PHASE TRANSITION SIGNATURES IN 
CLUSTER DISTRIBUTIONS 

A. Insensitive signatures 

In this section the insensitive features of the cluster 
distribution for each system are examined. It will be 
demonstrated that on this level of analysis each system 
exhibits behaviors that are consistent with systems which 
undergo a continuous phase transition. The conclusion is 
inescapable that this sort of analysis can yield necessary, 
but not sufficient, signals and no further insight to the 
mechanism behind multifragmentation. A deeper analy- 
sis will be necessary to distinguish those systems which 
undergo such a phase transition from those which do not. 



1. Fluctuations 

One of the most striking characteristics of systems un- 
dergoing continuous phase transitions is the occurrence 
of fluctuations that exist on all length scales in a small 
range of the control parameter. In fluid systems this 
was observed as critical opalescence, first noted by An- 
drews in the latter half of the 19 tft century |4S|| . Fluctu- 
ations in cluster size and the density of the system arise 
because of the disappearance of the latent heat at the 
critical point. This is illustrated in the FDM when the 
isothermal compressibility diverges at the critical point 
and small changes in pressure gives rise to great changes 
in the density. In the FDM as the volume and surface 
contribution to the free energy of cluster formation van- 
ishes the power law dominates and clusters of all length 
scales are observed M. 
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In a cluster distribution the most readily observed fluc- 
tuations are those in the size of the largest cluster. For 
each system the root mean square (RMS) fluctuations in 
the size of the largest cluster normalized to the size of the 
system, A (A max /Ao), have been calculated as a function 
of the system's control parameter. This measure of the 
fluctuations in the cluster distribution was first studied 
by Campi for gold multifragmcntation and percolation 
Q . Those results are replicated here for those two sys- 
tems. 

Figure la shows A (A max /Ao), as a function of pi for 
percolation. As expected for a system known to exhibit a 
continuous phase transition, the RMS fluctuations peak 
over a narrow range in the control parameter. The loca- 
tion of this peak provides a first estimate of the critical 
point; p c = 0.33 ± 0.01. See Table I. 




FIG. 1. Fluctuations in the normalized size of the largest 
cluster for: (a) percolation (Aq = 216) as a function of pi, 
(b) percolation as a function of m, (c) random partitions as 
a function of m and (d) Au + C multifragmentation as a 
function of m. Solid vertical lines indicate the guess for the 
critical point. 

Next the percolation lattice is examined using the mul- 
tiplicity of clusters, to, as an estimate of the control pa- 
rameter. This is done because in the case of nuclear mul- 
tifragmentation to is experimentally measurable. Figure 
lb shows much the same qualitative behavior as Figure 
la. The fluctuations peak over some narrow range of m 
and suggest the value of the multiplicity at the critical 
point, the critical multiplicity, to be m c = 55 ± 5. 

For random partitions a peaking behavior in the fluc- 



tuations of the size of the largest cluster as a function 
of m was observed, see Figure lc. These fluctuations 
can be understood as follows. At m = 1 there can be 
no fluctuations in the size of the largest cluster because 
of the dual constraints of event cluster multiplicity and 
the fixed number of constituents. As the multiplicity in- 
creases from one, the constraints ease and fluctuations 
in the size of the largest cluster grow. At the maximum 
possible multiplicity, i.e. when to is equal to the total 
number of constituents, the size of the largest cluster 
is constrained to be equal to unity. Thus, the fluctu- 
ations show a peak, but for a reason that has nothing 
to do with a continuous phase transition. Therefore it 
must be concluded that the observation of a maximum 
in the fluctuations in the size of the largest cluster is 
not sufficient to distinguish systems with and without 
critical behavior. On the other hand, the absence of a 
peak in fluctuations would indicate that the clusters of 
the system were not produced near a critical point. If 
the system's phase space has been fully explored, then 
the stronger statement that the system does not possess 
a critical point could be made. At this level of analysis 
the critical multiplicity of this system can be estimated 
to be to c = 5 ± 2. 

Finally, Figure Id shows the Au + C multifragmenta- 
tion data with the cluster distribution normalized to the 
size of the system, Ao(m). The fluctuations in the mass 
of the largest cluster exhibit a peak when plotted as a 
function of the event total charged particle multiplicity, 
m. This behavior is consistent with what is expected for 
a critical phenomenon. However, as illustrated above, 
it is far from conclusive. At this level of analysis the 
estimate for the critical multiplicity is m c — 31 ± 6. 

It is also possible to study the fluctuations in the aver- 
age size of a cluster. From the example of critical opales- 
cence it is clear that the greatest fluctuations in cluster 
size should occur at the critical point. To that end the 
quantity known as 72 is constructed again following the 
work of Campi |12| - fllTf . The variance in the mean 
cluster size, (A/), is defined as: 



(AfY 



(20) 



The average cluster size is given by the ratio of the first 
moment to the zeroth moment: 



Y j n Af A f /^2n Af =M 1 /M . 



(21) 



The first term in eq. ( f20| ) is just the ratio of the second 
moment to the zeroth moment. Therefore, the variance 
in the average cluster size can be written in terms of the 



>.th 



moments: 



a 2 = 



M2 
M 



( M > • 



(22) 



This quantity is directly related to Campi's 72 via: 



G 



TABLE I. Critical point determination 



Method / System 



Percolation (pi 



Percolation (m) 



Random Partitions 



An + C 



A (A max /A ) 
mi peak 

Fisher r-power law 
Scaling Function 
7-matching 



0.33±0.01 
0.28±0.03 
0.31±0.05 
0.34±0.03 
0.33±0.02 



55±5 
62±2 
57±3 
57±6 
49±1 



5±2 

5±2 

59±1 

10±2 

9±1 



31±6 
18±2 
22±1 
22±2 
21±2 



72 = 



{AY 



M 2 M 
Ml ' 



(23) 



which is easily measured and was coined by Campi as the 
reduced variance ]12|]. 

In a later paper, ]16[| , Campi discussed the differences 
in methods to measure 72. Specifically, the manner in 
which the fc^-moments are computed from the observed 
cluster distribution. One method is to measure the k th - 
moments on an event by event basis and then compute 
an average based on the control parameter, e.g.: 



1 



JV 



=1 



1 

N 



JV 



,(e)4)- ( 24 ) 



= 1 A f 



where N is the number of events at a control parameter 
value of e, and i denotes the i th event. This method of 
calculation of the fc^-moments will be termed averaging 
the sums and will yield: (72). 

The alternate method involves calculating an average 
cluster distribution at each value of the control parameter 
and then calculating the fc^-moments from the resulting 
average cluster distribution: 



1 N 

Mk = £ (n Af (e)) A) = B^E^ (e))A 



(25) 



This method of calculation will be termed summing the 
averages and will give: 7 2 . 

For quantities linear in tia, there is no difference in 
the two methods so that (Mf.(e)) = Mk(e). However, 
due to the dependence of 72 on the square of the first 
moment, there will be a difference in the two methods of 
calculation. Results for both methods for each system 
are shown in Figure 2. 

Of primary significance is the presence of a peak in 
both measurements of 72 for all systems. For an infinite 
system exhibiting critical phenomena, the location of the 
peak in 72 will coincide with the location of the critical 
point. For the percolation system Figures 2a and 2b show 
that both the location and magnitude of the peak in 72 
is dependent on the choice of calculation method. Solid 
lines indicate this measure of the critical point. For the 
random partitions, Figure 2c, the location of the peak 
in 72 shows no dependence on the method of calculation 
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FIG. 2. Fluctuations as measured by 72 for: (a) percola- 
tion as a function of pi, (b) percolation as a function of m, 
(c) random partitions as a function of m and (d) Au + C 
multifragmentation as a function of m (at low m 7 2 is af- 
fected by fission events not completely filtered). Open circles 
show (72), while filled circles show 7 2 (see text for details). 
Solid vertical lines indicate the guess for the critical point. A 
dotted horizontal line shows the value of 72 = 2. 



while the magnitude of the peak does. The gold multi- 
fragmentation data exhibit a dependence on the method 
of calculation both in the magnitude and location of a 
peak in 72. 

Having noted the peaking behavior of 72, the signif- 
icance of the amplitude of the peak is now addressed. 
It has been suggested that the height of the peak can 
be used to differentiate between the presence of a power 
law and that of an exponential: for a power law 72 > 2 
while for an exponential 72 < 2. This is not definitive 
proof of the existence of a continuous phase transition 
as other systems show power laws in the absence of such 
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a phase transition. All of the percolation figures show 
peaks above two, as do the multifragmentation data plots 
and the random partitions. However, the value of 72 de- 
pends on the size of the system in question [ fl6| |. For a 
percolation system with 64 sites, peaks in 72 under two 
are observed, see Figure 3a and 3c. Therefore, the crite- 
rion 72 > 2 is not sufficient to discriminate between those 
finite systems which do and those which do not posses a 
power law cluster distribution. 
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FIG. 3. Fluctuations as measured by 72 (see text for de- 
tails) for a percolation lattice of 64 sites: (a)open circles {72), 
filled circles 7 2 as a function of pi , (b)A72 as a function of pi, 

(c) open circles (72), filled circles 7 2 as a function of m and 

(d) A72 as a function of m. Solid vertical lines indicate the 
guess for the critical point. A dotted horizontal line shows 
the value of 72 =2. 

Finally the question of the difference between the al- 
ternative methods of calculating 72 is examined via: 
A72 = (72} — 72- It nas been suggested that a peak 
in the difference could indicate critical phenomena and 
the location of the critical point |l6) . Unfortunately, the 
cause of this peak is not well understood and vanishes 
at the limits of the system size: (0, 00). Figures 4a and 
4b do show peaks in A72 at some intermediate value of 
the control parameter for this percolation lattice of 216 
sites. However, as the size of the percolation lattice in- 
creases this signal vanishes [16[ . For a percolation lattice 
with 64 sites Figures 3c and 3d, respectively, look like a 
cross between the percolation (L — 6, m) results, Fig- 
ures 2b and 4b, and the random partition results shown 
in Figures 2c and 4c. This is believed to be due to the 



twin constraints of the multiplicity and the conservation 
of constituents imposed upon the system at the extremes 
in cluster multiplicity. Similar behavior is observed in 
the gold multifragmentation data in Figures 2d and 4d. 
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FIG. 4. The quantity A72 for: (a) percolation (A = 216) 
as a function of pi, (b) percolation as a function of m, (c) 
random partitions as a function of m and (d) Au + C multi- 
fragmentation as a function of m. Solid vertical lines indicate 
the guess for the critical point. 

Neither the 72 measure of fluctuations nor the observa- 
tion of fluctuations in the size of the largest cluster pro- 
vide definitive insight into the nature of the cluster pro- 
ducing mechanism. For both random partitions and per- 
colation 72 peaks at nearly the same value of the control 
parameter regardless of the method of averaging used. 
For the percolation system the value of pi at the peak in 
72 is close to the value of p\ where A(A max /Ao) is a max- 
imum. This coincidence does not hold for random par- 
titions; compare Figure lc and 4c. For both percolation 
(to) and multifragmentation, there is better agreement 
on the critical point from fluctuations and than from 72 
is computed via eq. (|24|). 



2. Divergences 

Another signature previously used to infer the exis- 
tence of a continuous phase transition from cluster dis- 
tributions is the observance of a peak in the second mo- 
ment [^lj, It has been pointed out that models 
with no phase transition can exhibit a peaking behavior 
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in the second moment H] . Figure 5 shows the behavior 
of the second moment for each of the systems examined 
in this work. In this figure, for the sake of illustration, 
the largest cluster has been excluded from the sum at all 
values of the control parameter. Each system shows a 
peak at some intermediate value of its control parameter. 
Table I lists the location of the second moment peaks. It 
is clear from the peak observed for the random partitions 
that it is possible to observe a peak in the second mo- 
ment for a non-critical cluster distribution. Thus this 
quantity cannot be used to distinguish between critical 
and non-critical systems. 
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FIG. 5. The second moment, M%, of the cluster distribu- 
tions for: (a) percolation as a function of pi, (b) percolation 
as a function of m, (c) random partitions as a function of m 
and (d) Au + C multifragmentation as a function of m. Solid 
vertical lines indicate the guess for the critical point. 

An issue with the use of the second moment's peaking 
behavior is the exclusion of the largest cluster from the 
sum in eq. (fill). Again, in the FDM formalism the sum 
runs over all clusters in the gas. On the liquid side of the 
critical point a gas exists in addition to a liquid drop. 
Thus, the largest cluster represents the bulk liquid. On 
the gas side of the critical point there is no liquid drop 
and the largest cluster is merely the largest gas particle. 
With this understanding it is clear that the largest clus- 
ter should be omitted from the summation in the second 
moment only in the liquid region, whereas the summa- 
tion should run over all clusters in the gas region. For a 
proper construction of the second moment, knowledge of 
the location of the critical point is required. In the ther- 



modynamic limit of infinite system size, exclusion of the 
largest cluster makes little difference. However in small 
systems the proper construction of the second moment is 
crucial if critical behavior is to be observed in ref. IpOl . 



See ref. 61 for an example of the improper construction 
of the second moment. 



3. Campi plots 

Plots of the natural log of the normalized size of the 
largest cluster, \n(A max j Aq) , versus the natural log of the 
second moment, ln(M2), were first presented by Campi 
in a comparison of gold multifragmentation and percola- 
tion [[fl]]. Figure 6 shows the resulting plots for each of 
the systems discussed in this paper. In each plot there is 
a liquid leg for the largest A max and small M2 and a gas 
leg for smaller A max and mid-range values of Mi . That 
similar behavior is observed for all systems is in indica- 
tion that this is a necessary, but not sufficient, signal for 
critical behavior. 
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FIG. 6. Campi plots of the cluster distributions for: (a) 
percolation as a function of pi, (b) percolation as a function 
of m, (c) random partitions as a function of m and (d) Au + 
C multifragmentation as a function of m. 



4- Rise and fall of intermediate mass fragments 

In many nuclear multifragmentation studies the term 
intermediate mass fragment, IMF, has been defined as a 
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cluster which has a charge between 3 < Zf < 30. For 
the percolation system presented here a charge has been 
assigned to each cluster by multiplying the number of 
constituents in the cluster by the charge to mass ratio of 
a gold nucleus. For the random partitions the number of 
constituents is used as the charge. Since the definition of 
an IMF is arbitrary it makes little qualitative difference 
what range in some measure of the cluster size is used. 

Aside from the equilibrium arguments made by the 
ALIDIN group Q- Q, little insight towards the pres- 
ence or absence of a continuous phase transition is gained 
from a simple plot of the average number of IMF's, 
(Mi m f) versus the control parameter. Figure 7 shows the 
results for the systems discussed in this work. Each sys- 
tem shows a peak in (M im f) at some intermediate value 
of the control parameter. Comparing the peak position 
in Figure 7 to the values listed in Table I shows that there 
is little correspondence between the numerous proposed 
methods for locating the critical point. The arbitrary 
nature of the definition of an IMF makes it unlikely that 
the peak in (Mi m f) occurs at the critical point. To some 
degree the rise and fall feature is due to the constraint 
of a fixed number of constituents. It it obvious that a 
the extreme values of the control parameter, the num- 
ber of IMFs must diminish, while at intermediate values, 
it must be at least as great. Thus, the occurrence of a 
peak at some intermediate value of the control parameter 
is expected. 
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FIG. 7. Average number of intermediate mass clusters, 
(Mimf), for: (a) percolation as a function of pi, (b) percola- 
tion as a function of m, (c) random partitions as a function 
of m and (d) Au + C multifragmentation as a function of m. 



5. r e ff -minimum 

With the first observation of a power law in the nuclear 
multifragmentation yield distribution |Q, |7| it became 
a common analysis tool to fit cluster distributions to a 
power laws and extract exponent values. In an effort to 
make this a more quantitative analysis the value of the 
extracted exponent, r ( . //. was examined as a function of 
some control parameter that was experimentally or nu- 
merically accessible. It was assumed that at the critical 
point the value of r e ff should attain a lower value than 
fits which were performed away from the critical point 
p2| - |p7fl . The logic of this assumption was based upon 
the idea that at low temperatures a system has few small 
clusters, so the power law should be steep, leading to a 
high r e ff value. At high temperatures there are many 
small clusters and little else, which is reflected in a high 
value of T e ff and a steep power law. At the critical point 
clusters on all length scales appear and the power law is 
shallow with a lower value of T e //. In this analysis the 
largest cluster was generally omitted from the fitting pro- 
cedure and both the constant of proportionality and r e / / 
were allowed to vary independently. Many investigations 
of nuclear multifragmentation, both theoretical and ex- 
perimental, employed this method of analysis p2[ - |67t). 



There are two flaws in this analysis method. The first 
is the use of a two parameter fit for the power law. Al- 
lowing both the overall normalization of the power law 
and the exponent to vary independently is in conflict to 
the scaling assumptions underlying the FDM as shown in 
eq.'s (pp and <^). A proper fit for a power law within the 
context of the FDM should be based on single parameter. 
Furthermore, the cluster distribution must be normalized 
to the size of the system as was outlined in III. With- 
out this normalization, which requires knowledge of the 
system's size, power law fits lose much of their ability to 
contribute useful information to the presence of critical 
phenomena. 

Leaving aside for a moment that the execution of the 
r e //-minimum analysis violates the scaling assumptions 
of the FDM, the signal of a minimum in the cluster yield 
power law will be examined. A two parameter fit for r e f f 
searches for the minimum in an effective exponent which 
is defined as floSll - [frol: 



T eff = - 



d In n J 4 / (e) 



d\nA 



f 



(26) 



If it is assumed that the system under study follows a 
power law in the cluster yield at the critical point, and 
away from the critical point the cluster yield is affected 
by a scaling function such as in eq. (Q), then: 
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= T-A 



d In/ 



(27) 



The minimum in r e f / can be found by differentiating eq. 



dT e ff d dhif 

= -Afj^^r^- = 0. 



de 



de OA 



f 



(28) 



This indicates that the location of the minimum in r e f f 
is dependent on the form of the scaling function, /. As- 
suming the scaling function has the form of eq.(^|) then 
the minimum in r e ff will be at e = B/2A'j , and not at 
the critical point e c = 0. 

Despite the flaws in the r e ff -minimum analysis it 
is of interest to examine the results for the systems 
discussed in this paper. Figures 8 through 11 show 
the results for a two parameter fit to the cluster dis- 
tribution for percolation (probability and multiplicity), 
random partitions and gold multifragmcntation, respec- 
tively. For all systems, the cluster distributions were fit 
at each value of the control parameter. Only clusters 
with -0.02 < (Af/Ao) < -0.22 were included in the fits. 
The first three systems weighted xl with errors associ- 
ated with n,A f (e) while the xl for the gold multifragmen- 
tation cluster distributions were weighted with errors on 
both nA f {e) and Af. 

For the percolation (pi) a minimum in r e ff was ob- 
served at pi = 0.3 with xl = 2.3, q = 0.214 ± 0.005 
and t = 2.19 ± 0.01; shown in Figure 8a, b and c by 
the dotted lines. However, at p\ = 0.33 the xl = 1-02 , 
q = 0.181 ± 0.003 and r = 2.27 ± 0.01; shown in Figure 
8a, b and c with the dashed lines. Based on a goodness 
of fit comparison, the latter value of pi is a better choice 
for the critical point as the cluster distribution is better 
fit by a power law. This result is in agreement with the 
analytic discussion of r e ff above, namely that a mini- 
mum in r e ff is a poor indicator of the critical point. If 
the results for pi = 0.33 are compared to the center of the 
r e ff, pi — 0.28, the differences in the xl an d go results 
are even more striking. 

Similar results were seen for percolation (m), see Fig- 
ure 9. Here the minimum in the r e y/-well yielded worse 
results for both xl an d 1o than does the choice of the 
critical point based on a choice from the xl ~ 1 region 
where there is good agreement between the fitted go and 
the value computed using eq. ([)]) and the canonical r 
value for three dimensional percolation. 

Significant differences between percolation and random 
partitions are observed in this analysis as seen in Figure 
10. The solid lines show the r e ff and go values for sys- 
tems in the three dimensional Ising universality class, 
while the dashed line shows the r e f f and qo for three di- 
mensional percolation. The first noticeable difference is 
a lack of a valley shape in the plot of r e f f versus control 
parameter, see Figure 10b. The value of T e // is below 
2.2 for all but m > 60. Next is the lack of a region in 
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FIG. 8. Results for r e //-minimum analysis on percolation 
as a function of pi: (a) xl-> (b) fitted r e // and (c) qo- The 
vertical dotted line show the location of the minimum in r e //. 
The vertical dashed line shows one instance of a better fit 
based on xl- The horizontal solid lines show the accepted 
values of r and go for percolation. 



m where xl ~ 1 (other than at m = 2), Figure 10a. All 
fits yield large xl values indicating poor fits to the clus- 
ter distribution by a power law for the range of clusters 
examined. 

The gold multifragmcntation data show results similar 
to those of percolation (m) . Here the cluster size is mea- 
sured in terms of the nucleon number and the cluster 
distribution is normalized to the mass of the gold pro- 
jectile remnant. Figure lib shows a valley in r e // as a 
function of m, albeit one with a shallow and question- 
able upwards slope at high m. Figure 11c shows fitted 
values of go that coincide with canonical values. Figure 
11a shows a region of low xl values followed by steadily 
increasing values. If no knowledge of the go and r values 
is assumed, then this analysis shows no definitive signals. 
The r e //-valley shows a broad minimum in xl thus no 
one value of m can be selected for the critical point based 
on goodness of fit arguments. At best one could argue 
for the neighborhood of the critical point and a value of 
go and r in some broad range. 
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FIG. 9. Results for r e //-minimum analysis on percolation 
as a function of m: (a) xt, (b) fitted r e // and (c) go- The 
vertical dotted line show the location of the minimum in r e f / • 
The vertical dashed line shows one instance of of a better fit 
based on xt- The horizontal solid lines show the accepted 
values of r and go for percolation. 



6. Conclusion 

The analysis presented in this section shows that many 
proposed indicators of critical behavior are inconclusive. 
All of the considered systems show similar signals which 
are qualitative in nature and open to interpretation. It 
is therefore impossible, based solely on this level of anal- 
ysis, to make a definitive conclusion as to the presence 
of a continuous phase transition in any of these systems. 
What is needed is an analysis or set of analyses that more 
clearly differentiates between systems with and without 
critical behavior. 



B. Sensitive signatures 

1. The Fisher t -power law and the critical point 

In this section the cluster yields are fit to a power law 
in a manner consistent with the FDM formalism. As 
with the two parameter fits the yields for clusters with 
~0.02 < (Af/A ) < ~0.22 were fit at each value of the 
control parameter. However, only a single parameter, the 
value of r, was allowed to vary to minimize the xt of the 
fit. The value of the normalization, qo, was determined 



FIG. 10. Results for r e //-minimum analysis on random 
partitions: (a) xt, (b) fitted r e ff and (c) qo- The horizontal 
solid (dashed) lines show the accepted values of r and qo for 
percolation (3D Ising). 

via the Riemann ^-function in eq. (j^). As suggested by 
Fisher ||, the value of r was constrained to be between 
2 and 3 so that the sum in the ^-function converges. 

If the cluster distribution is well described by the FDM, 
then at the critical point the fit to a single parameter 
power law should show a minimum in xt- Away from 
the critical point the power law is modified by a scal- 
ing function with a form similar to that given in eq. 
Therefore, fits to a single parameter power law should be- 
come worse as the modification from the scaling function 
increases away from the critical point. 

Figure 12 shows the results for the percolation sys- 
tem with pi as the control parameter. In Figure 12a, a 
minimum in xt i s observed for fits in the mid pi range. 
This minimum indicates the location in pi of the cluster 
yield distribution which is best fit by a single parame- 
ter power law as suggested by the FDM formalism. By 
this estimation the critical point for this 216 site perco- 
lation lattice is at p c — 0.31 ± 0.05 with r = 2.2 ± 0.1, 
<7o = 0.20 ± 0.01 and a xt = 1-62. The canonical values 
of r and qo are not extracted due in part to unavoidable 
finite size effects, and to the binning of cluster yields to- 
gether over a range of 0.01 in pi, which causes the true 
cluster distribution at the critical point to be contami- 
nated by distributions at other values of the control pa- 
rameter. In spite of these difficulties, the signature of 
the critical behavior suggested by the FDM formalism 



12 



10 - 



(a) 




10 E 



1 





• • * • 

• • • 
• • 

(a) 









3.2 

3 
2.8 
2.6 
2.4 
2.2 

2 
1.8 



0.1 



0.2 



0.3 



0.4 



0.5 



0.6 



0.7 



h 


(b) 




\ 


- 


l i i i l i i i 


0.2 


0.4 


0.6 0. 






Pi 



0.8 

Pi 



« 


(c) 




f 
* 
* 
• 


* -Ilffllfl 

, 1 , , , 1 , , , 1 



0.2 0.4 0.6 0.8 
Pi 



FIG. 11. Results for r e //-minimum analysis on Au + C 
multifragmentation: (a) xt, (b) fitted r e // and (c) go- The 
horizontal solid (dashed) lines show the accepted values of r 
and qo for percolation (3D Ising). 



FIG. 12. Results for the full Fisher r-power law analysis 
on percolation as a function of pi : (a) \t , (b) fitted r and (c) 
go . The vertical dashed line shows the estimate of the critical 
point based on a best fit based on xt- The horizontal solid 
lines show the accepted values of r and go for percolation. 



is unmistakable. The location of the critical point de- 
termined here is consistent, at the 10% level, with the 
insensitive signatures presented in the previous section. 
See Table I. Figure 16a shows the best fit power law. For 
the percolation system clusters consisting of a single site 
are excluded from the fitting procedure. It is accepted 
that those clusters reflect the effects of the finite size of 
the system to a higher degree than larger clusters. Clus- 
ters with Af < 53 were included in the fit. The largest 
cluster from each event was excluded from consideration 
when generating the average cluster distribution in keep- 
ing with the FDM formalism. Figure 16a shows the data 
for the entire cluster distribution in open circles. It is 
clear from this figure that the majority of the cluster dis- 
tribution was used in the power law fit and further, that 
the exclusion or inclusion of the larger clusters has al- 
most no effect on the results of this procedure. In short, 
the extracted parameters, namely r, qo and p c , do not 
depend on the fit range. 

Figure 13 shows the results of the single parameter 
fit analysis when applied to the same percolation system 
but with the cluster multiplicity used as a measure of the 
control parameter. Again there is a minimum in the xt 
values at some intermediate value of the control parame- 
ter which indicates that m c — 57 ± 3 with r = 2.2 ± 0.1, 
qo = 0.20 ± 0.01 and a, xl = 0.72. Note the consistency 
between these values of q and r and those obtained with 



pi following this method. The location of the critical 
point determined here is consistent, at the 10% level, 
with the insensitive signatures presented in the previous 
section. The lower xt value is due to the finer bins over 
which the cluster distributions were grouped. Figure 16b 
shows the best fit power law. Here only clusters of size 
A = 1 and size A = A max were excluded from the fitting 
procedure. 

From Figures 12 and 13 it could be argued, based on 
the best agreement between the fitted r and the accepted 
three dimensional percolation value, that there are bet- 
ter choices for the critical point than those quoted above. 
However, those arguments assume knowledge of the value 
of t as an input. The use of the location of the best fit 
to a single parameter power law as an indicator of the 
critical point makes no assumption regarding the value 
of t and is a test of the FDM formalism in which only 
the range of r is suggested: 2 < r < 3. The values 
of t and go are outputs rather than inputs of this anal- 
ysis. Much of the following analysis presented in this 
paper follows the same philosophy. That is, the analysis 
is designed to test the cluster distribution in question for 
behavior consistent with the FDM formalism. The val- 
ues of quantities, such as critical exponents, are results 
of the analysis method and are in no way selected for 
on the basis of their particular values. Agreement be- 
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FIG. 13. Results for the full Fisher r-power law analysis 
on percolation as a function of m: (a) \ v , (b) fitted r and (c) 
go • The vertical dashed line shows the estimate of the critical 
point based on a best fit based on xl- The horizontal solid 
lines show the accepted values of r and go for percolation. 



tween exponent values determined by this procedure and 
the canonical values in various universality classes is then 
significant because the values of the exponents are deter- 
mined solely by the behavior of the cluster distributions 
so analyzed. 

The results of the single parameter power law fits for 
the random partitions are presented in Figure 14. There 
is a minimum in the xl value at to = 59. However 
xl = 10.83, which is an order of magnitude above the 
percolation results, should not be used as an indication 
of a good fit of the cluster distribution by a single param- 
eter power law. The location of the xl minimum is also 
in disagreement with the insensitive signatures presented 
in the last section. Here only clusters of size Af = 1 and 
size Af — A max were excluded from the fitting proce- 
dure. 

Figure 15 shows the results of this analysis applied to 
the gold multifragmentation data. As with the percola- 
tion results, the xl shows a minimum that drops nearly 
two orders of magnitude from the peaks for high and low 
to to the valley at a mid range value of to, see Figure 15a. 
In the context of the FDM analysis this result suggests 
that the critical point is located at m c = 22 ± 1 with 
t = 2.2 ± 0.1 and q a = 0.18 ± 0.01 and xl = 2-70. The 
best fit power law is show in Figure 16d. An uncertainty 
of one unit of multiplicity is assigned to m c to take into 



FIG. 14. Results for the full Fisher r-power law analysis 
on random partitions: (a) xt, (b) fitted r and (c) go- The 
horizontal solid (dashed) lines show the accepted values of r 
and qo for percolation (3D Ising). 



account the relatively low values of xl of the neighboring 
fits. 

For the above fits to the gold multifragmentation data 
the xl is weighted by the errors in both ua s and A / . The 
fitting procedure has also been performed with no error 
weighting on xl and with errors only in ha s for weight- 
ing. Both analyses shows results that were not signifi- 
cantly different from those quoted here. As mentioned 
previously, clusters with Zj = 2 are created in both the 
prompt first stage and in the multifragmcnting of the gold 
nuclear remnant. The prompt Zf = 2 clusters have been 
excluded from the filtered gold multifragmentation anal- 
ysis. However, as a further test of the single parameter 
power law fit, only clusters with 3 < Zf < 16, i.e. clus- 
ters with no contamination from the prompt first stage, 
were included in a repeat of this analysis. Again the re- 
sults show practically the same behavior as those shown 
here. As yet another test, clusters with 2 < Zf < Z max 
were included in the fitting procedure, and again the re- 
sults showed no difference from those presented here. Fi- 
nally clusters with 3 < Zf < Z max were included in the 
fitting procedure, and again the results showed no differ- 
ence from those presented here. The data always showed 
a deep valley in the \l versus to plot which indicated that 
the location of the critical point was m c ~ 22 and that 
t ~ 2.2, q <~ 0.18 and 1 < xl < 4. This analysis shows 
that the value of r and the location of m r are not sensi- 
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FIG. 15. Results for the full Fisher r-power law analysis 
on Au + C multifragmentation: (a) xi> 3 ) fitted r and (c) 
qo. The vertical dashed line shows the estimate of the crit- 
ical point based on a best fit based on xt- The horizontal 
solid (dashed) lines show the accepted values of r and qo for 
percolation (3D Ising). 
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FIG. 16. Best fit Fisher r-power laws for: (a) percolation 
(pi), (b) percolation (m), (c) random partitions and (d) Au 
+ C multifragmentation. Solid circles were included in the 
fits, open circles were excluded. Solid lines show the best fit 
power law, dotted lines indicate error in the power law based 
on errors in the fitted parameters. 



tive to the fit region. The behavior of the data show this 
clearly. See open circles in Figure 16d. 

The single parameter power law analysis of the cluster 
distributions of the various systems produced the first 
result which can differentiate between systems that fol- 
low the FDM formalism and systems that do not. The 
differences between Figures 12a, 13a, 15a and Figure 14a 
are clear and distinct. For both percolation and gold 
multifragmentation the behavior of xt is exactly what is 
predicted by the FDM formalism for continuous phase 
transitions. Far from the critical point the cluster distri- 
bution is fit poorly by a single parameter power law due 
to the influence of a scaling function where volume and 
surface effects overwhelm the underlying power law. At 
the critical point, where the influence of the scaling func- 
tion vanishes, the cluster distributions are well described 
by a single parameter power law with an exponent value, 
r ~ 2.2 and thus qo ~ 0.2, in keeping with what is ex- 
pected for many universality classes. This fitting pro- 
cedure does not merely search out a cluster distribution 
which is well fit by a power law, but finds the cluster dis- 
tribution which is well fit by the FDM formalism. This 
is achieved via the coupling between the exponent r and 
the normalization factor qq. See eq. (Q). The random 
partitions fail to produce such signals. This is expected 



as that system does not obey the FDM formalism and 
thus should not show the same signals as systems that 
are known to follow the FDM such as percolation. This 
analysis of the cluster yield of gold multifragmentation 
yields a signal that is suggestive of critical phenomena. 

2. The critical exponent a 

In section III it was shown that in the context of the 
FDM the surface of a cluster makes a contribution to the 
free energy of cluster formation, via the scaling function 
f(z), that depends on the number of constituents of the 
cluster raised to the power a. See eqn's ([!]), and 
The behavior of the order parameter suggests that the 
scaling function f{z) has a maximum [ p2[ . At the maxi- 
mum of the scaling function, f max (zmax)} the production 
of Af sized clusters is greatest: 

n%; x (e max ) = qoAj T f(z max ). (29) 

The argument of f m ax is: 

Zmax — AfCrnax, (*^J 

where the value of z max depends on the specific details of 
the system in question [p3|. Rearranging eq. (|30|) yields: 
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^max — Z max A j . 

Thus tmax-, the value of the control parameter at which 
the greatest number of clusters of size Af are produced, 
is related to the cluster size through a simple power law 
with exponent a. The exponent a can then be deter- 
mined from knowledge of the location of the critical point 
and the value of the control parameter at the greatest 
production of clusters of size Af. 

The location of the critical point was determined in the 
search for the Fisher r-power law in section IV-B-1 and 
will be used here to determine the a. The value of the 
control parameter which yields the greatest production 
of each Af cluster size was determined from the peak 
location in a plot of riA f (e) versus the system's control 
parameter. See Figure 17. 




FIG. 17. Rise and fall in the production of individual clus- 
ters for: (a) percolation as a function of pi , (b) percolation as 
a function of m, (c) random partitions as a function of m and 
(d) Au + C multifragmentation as a function of m. Symbols 
(filled O) D i A) v) indicate clusters of size (charge) 3, 5, 7 
and 10 for percolation and random partitions (multifragmen- 
tation). 

For each system at each cluster size plots such as those 
shown in Figure 17 were used to determine the location 
of the peak of n^ f (e). For example, in percolation (pi), 

trnax = (Pc ~ Pmax) / Pc, of points {n Af {f),pi) for 

a particular Af were fed into a SPLINE routine J72[ . 
Input pairs were then smeared by assigning SriA f (e) as 
the standard deviation of a gaussian centered on riA^e). 
Output of the SPLINE routine was used to interpolate 



the behavior of a smooth curve between the pairs of input 
points. Stepping along the interpolations in increments 
much smaller than the separation of the input pi, a max- 
imum of riAj (e) was determined and p m ax was recorded. 
This process was repeated thousands of times for each 
cluster size and lead to an estimate of p m ax ± Sp max as a 
function of Af. 

Using p max ( Af ) ± Sp max {Af ) and the value of p c ± Sp c , 
from the Fisher r-power law determination process, the 
value of the exponent a was determined by taking the 
slope of hi(e m ax) versus \n(Af). The value of z max was 
determined by exponentiating the offset. The value of p c 
was varied uniformly throughout the range suggested by 
Sp c and tens of fits were made with varying starting and 
ending points in Af of the fitting region. The final value 
of o i 5<T and z max ± Sz max are the average and RMS 
values resulting from all the fits. 

Results of this analysis performed on percolation (pi) 
are shown in Figure 18a. Here the value of the con- 
trol parameter that coincides with the maximum in pro- 
duction of clusters of size Af, e max , is plotted against 
the cluster size. Results of the average power law fits 
to eq. ( |3T| ) are plotted as a solid line. The agreement 
between the values returned by the procedure discussed 
above, a = 0.52±0.02 and z max = -0.89±0.03, and 
the accepted values for three dimensional percolation, 
a = 0.45 and z max = —0.8 Q, establishes the reliability 
of this exponent extraction method. The analysis differs 
in method from previous efforts on percolation lattices 
0, Q but not in result. 

The next test of this analysis is to extract a value of a 
from percolation (to). In order for this procedure to be 
useful on multifragmentation data it must be shown that 
the exponent a can be determined using cluster multiplic- 
ity as the control parameter. To that end the multiplicity 
at which the production of each cluster size is maximal, 
fnmax was determined via the procedure described previ- 
ously. Using the value of m c ±Sm c determined via search- 
ing for the Fisher r-power law and m max the exponent 
a was determined by taking the slope of ln(e maa; ) versus 
In(Af). The value of m c was varied uniformly through- 
out the range suggested by 5m c and several plots were 
made with varying starting and ending points in Af of 
the fitting region. The value of z max was determined by 
exponentiating the offset. Results of the average power 
law fits to eq. ([u]) are plotted as a solid line in Figure 
18b. The agreement between the values returned by this 
procedure, a = 0.52±0.02, the value for a quoted in the 
above paragraph and the accepted values for three di- 
mensional percolation again establishes the reliability of 
this exponent extraction method and shows that the use 
of to as a control parameter is acceptable. 

The value of z max = —2.4 ±0.1 extracted for perco- 
lation as a function of multiplicity is different from the 
value quoted above, z max = — 0.89±0.03, for the per- 
colation system as a function of probability. This is a 
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FIG. 18. Results of the a extraction procedure for: (a) 
percolation as a function of pi, (b) percolation as a function 
of m, (c) random partitions as a function of m and (d) Au + 
C multifragmentation as a function of m. 

result of changing the measure of the control parameter 
from probability to multiplicity. This difference was ob- 
served in previous percolation efforts jnj and explained 
therein. A plot of e(pi) against e(m) show that z max (pi) 
and z max (m) map to each other. See Figure 9 of ref. Rnj . 

Clusters from the random distribution were also sub- 
jected to this analysis. Due to the failure of the search 
for the Fisher r-power law in the random partitions, the 
value of m c determined in the analysis of the gold multi- 
fragmentation was used, m c — 22 ± 1. The value of the 
cluster multiplicity for maximum production of Af sized 
clusters was determined in the same manner as with the 
percolation system. The flatness of the riA s (e) versus 
to curves, see Figure 17c, makes finding a unique value 
of m max impossible. This is reflected in the large error 
bars on e max ± 5e max seen when plotted against Af in 
Figure 18c. The value of m max reported by the peak find- 
ing procedure employed here reflects, approximately, the 
mid point of the multiplicity range of ua, (e) for a par- 
ticular A f . Coupling the to c from the filtered gold mul- 
tifragmentation data with the m max and fitting hi(e max ) 
versus \n(Af) gave a = 0.4 ± 0.2 and z max = —2.0 ± 0.8. 
However, it is clear when comparing the resulting average 
fit for the random partitions shown in Figure 18c with ei- 
ther of the percolation results shown in Figure 18a and 
18b that the a resulting for the random partitions clus- 
ter distribution is meaningless. This is to be expected as 



the framework of the FDM, used in the extraction of the 
exponent er, is meaningful only when applied to systems 
which undergo a continuous phase transition. The fail- 
ure of this analysis on this system is expected based on 
the basis of the failure of the analysis in the preceding 
section that aimed to find the Fisher r-power law and 
the critical point. 

Results for the extraction of a from the gold multi- 
fragmentation data have been published previously [||, 
pOf . In those analyses the largest cluster was excluded 
from consideration at every value of the control param- 
eter. This is at odds with the formalism of the FDM 
where the sums excludes the largest cluster for e > and 
include the largest cluster for e < 0. 

The previous analyses yielded values of a = 0.68 ±0.05 
and 0.65 ± 0.06 for work with the un-normalized charge 
distribution and normalized mass distribution respec- 
tively. When this analysis was redone using formalism 
of the FDM, i.e. the largest cluster excluded on one side 
of the critical point (liquid) and included on the other 
side (gas), the values of a were reduced by approximately 
50%: a — 0.32 ± 0.02. In the case of percolation the dif- 
ference introduced in the value of a when following the 
FDM formalism (as was done above) or not (as was the 
case in ref. (7^]) is on the order of a few percent. This 
is the first notable difference observed in the qualitative 
behaviors of percolation cluster distributions and gold 
multifragmentation cluster distributions. 

One source of this differing behavior is the changing 
mass of the system. For gold multifragmentation, from 
A - 194 at low m to A ~ 92 at high to Q, while the 
system size is constant for percolation. For gold multi- 
fragmentation effects of the finite size of the system are 
felt more at high multiplicities than low. Since the per- 
colation system size is constant, finite size effects are felt 
more evenly. 

It is the higher values of to where cluster production 
peaks in multifragmentation. The finite size of the sys- 
tem limits the size to which a cluster can grow. Thus the 
number of clusters of size At, function of to is 

contaminated when the largest cluster, A max , is included 
in a plot of ua s versus Af because A max would like to 
be larger, but finite size effects limit the size A max can 
attain. Therefore, one method to account for this effect 
is to exclude A max from the cluster distribution at large 
to values where these effects are largest This was done 
for the gold multifragmentation data. 

The multiplicity at which the production of each clus- 
ter size is maximal, m max was determined via the proce- 
dure described previously. The value of m c determined in 
the Fisher r-power law analysis was used, m c = 22 ± 1 . 
The value of m c was varied uniformly throughout the 
range suggested by 5m c and several plots were made with 
varying starting and ending points in Af of the fitting re- 
gion. The exponent a was determined by taking the slope 
of hi.(e max ) versus ln(A/) and the value of z max was de- 
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termined by exponentiating the offset. The results were 
a = 0.64±0.05 and z max = -6.0±0.8, see Table II and 
III. The average power law fits are shown in Figure 18d. 



3. The scaling function f(z) 

With the critical point (p c or m c ), r, qg and a deter- 
mined and assuming coexistence, g = 1, it is possible to 
find the scaling function by rewriting eq. (o) as 



n Af (e)/q A f T = f(z). 



(32) 



Doing this has the effect of appropriately scaling tla s (e) 
and collapsing the data onto a single curve. Figure 19 
shows the results of this sort of scaling. 

In Figure 19a, b and d, the data from percolation (pi 
and m) and multifragmentation, respectively, show col- 
lapse onto a single curve for a wide range in cluster size 
and over nearly the full range in control parameter. Ran- 
dom partitions shows no such collapse, see Figure 19c. 

As a demonstration of this type of scaling the same 
data has been scaled in the same fashion, but with a 
different choice of the critical point. Figure 20 shows 
the systems using a critical point with a value of half of 
the critical point determined via the Fisher r-power law, 
while Figure 21 shows the same analysis with a value of 
twice the critical point determined via the Fisher r-power 
law. A visual inspection of Figures 19, 20 and 21 reveals 
the greatest data collapse occurs when the choice of the 
Fisher r-power law critical point is used, at least for the 
percolation (pi and m) and multifragmentation systems. 
Random partitions show no such collapse. Using differ- 
ent values of r and a in this scaling analysis of random 
partitions does not significantly alter the data collapse. 
In Figures 19, 20 and 21 error bars on the data points 
are not shown for the sake of clarity. The size of the er- 
ror bars reflect the scatter of the data and are larger for 
larger negative values of z since there were lower statis- 
tics for higher multiplicity events |51|. 

Figure 22 shows a quantitative measure of the data col- 
lapse from this scaling analysis. For a number of different 
choices of control parameter scaling plots, as in Figures 
19 through 21, were made. Each plot was binned along 
the abscissa and the RMS fluctuations for each bin were 
calculated. The RMS fluctuations in all bins were then 
summed and plotted as a function of the choice of criti- 
cal point. See Figure 22. In the percolation (j>i and m) 
and multifragmentation systems the data shows the most 
collapse in the neighborhood of the Fisher r-power law 
critical point. No such behavior is observed in the ran- 
dom partition system. This analysis serves as another, 
albeit crude, estimate of the location of the critical point. 
Table I lists the results. 

The scaled data were used to determine the functional 
form of the scaling function by fitting the data with an 
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FIG. 19. Results of the scaling function analysis with the 
critical point equal to the value as determined from the Fisher 
r-power law analysis for: (a) percolation as a function of z(pi), 
(b) percolation as a function of z(m), (c) random partitions 
as a function of z(m) and (d) Au + C multifragmentation 
as a function of z(m). Solid curves show the fitted scaling 
function. Smaller sized points show smaller sized clusters. 
Cluster sizes are most evident in (d) where at the bottom 
in red are Aj = 2 sized clusters and at the top in orange 
are Af = 14. For percolation the colors cover a proportional 
range in cluster size while for multifragmentation the range 
is: red Zf — 3 to orange Zf — 14. 



empirical parameterization consisting of two gaussians 
instead of the single gaussian in eq (Q): 

/(z) = ai exp[-i(i-^) 2 ] + 
2 c\ 

a 2 eM-\C—) 2 ]- (33) 
2 c 2 

This was suggested by the asymmetry of the percolation 
(jpi) data, Figure 19a, and is consistent with a simplified 
version of corrections to scaling |6§| ] as discussed in sec- 
tion V. Figure 19 shows the resulting fits for all systems. 
Fit parameter values can be found in Table IV. Errors on 
the parameters of the fits, e.g. a\ etc., reflect the change 
in those parameters when the range of clusters included 
was changed, e.g. clusters with Zf — 2 were included 
or excluded and so on, and the weighting on the fit was 
changed, e.g. xt is unweighted, weighted with errors on 
nA f {e)/qoAj T or with errors on riA f (e)/qoAj T and e. 
The scaling function for percolation (pi and m) deter- 
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TABLE II. Critical exponents 



Exponent / System 3D Percolation Percolation (pi) Percolation (m) Random Partitions Au + C 3D Ising 



T 


2.18 


2.2±0.1 


2.2±0.1 


2.0±0.1 


2.2±0.1 


2.21 


a 


0.45 


0.52±0.02 


0.52±0.02 


0.4±0.2 


0.64±0.05 


0.64 




1.82 


1.5±0.2 


1.5±0.2 


2.5±1.2 


1.3±0.2 


1.23 


7+ (matching) 




1.8±0.2 


1.64±0.04 


0.4±0.1 


1.4±0.3 




7_ (matching) 




1.8±0.2 


1.7±0.1 


0.5±0.2 


1.4±0.3 




(7) (matching) 




1.8±0.2 


1.67±0.05 


0.5±0.1 


1.4±0.3 




A7 (matching) 




0.0±0.3 


0.06±0.1 


0.1±0.2 


0.0±0.4 




v (hyperscaling) 


0.87 


0.77±0.07 


0.77±0.07 




0.63±0.07 


0.63 



TABLE III. Critical amplitudes 



Amplitude / System 


Percolation (pi) 


Percolation (m) 


Random Partitions 


Au + C 


go 


0.20±0.01 


0.20±0.01 


0.07±0.01 


0.18±0.01 


Zmax 


-0.89±0.03 


-2.4±0.1 


-2.0±0.8 


-11.0±2.0 


Zmax (scaling fen) 


-0.9±0.1 


-2.5±0.5 


-2.0±8.0 


-9.0±2.0 


T + (scaling) 


0.9±0.1 


4.9±0.5 


3.5±0.5 


19.0T3.0 


T_ (scaling) 


0.17±0.05 


0.3±0.1 


1.0±0.2 


0.24±0.05 


r+/r_ (scaling) 


5.0±2.0 


16.0±6.0 


3.5±0.9 


80.0±20.0 


r+ (matching) 


1.0±0.3 


7.0±0.6 


26.0±2.0 


26.0±9.0 


T_ (matching) 


0.08±0.07 


0.28±0.04 


3.7±0.4 


0.27±0.06 


r+/r_ (matching) 


13.0±12.0 


25.0±4.0 


7.0±1.0 


100.0±40.0 


r+ (C.T.S. 3dP) 


0.9±0.1 


6.4±0.5 


170.0±20.0 


30.0±5.0 


r_ (C.T.S. 3dP) 


0.06±0.01 


0.38±0.04 


0.08±0.01 


0.08±0.01 


r+/r_ (c.t.s. 3dP) 


15.0±3.0 


17.0±2.0 


2100.0±400.0 


380.0T60.0 


r+ (C.T.S. 3dl) 






140.0±10.0 


55.0±5.0 


r_ (C.T.S. 3dP) 






0.1±0.01 


0.28±0.05 


r+/r_ (c.t.s. 3di) 






1400.0±200.0 


200.0±40.0 


j^scaled 








23.0±3.0 


j^scaled 








0.14±0.05 


j^scaled ^-pscaled 








170.0±60.0 



TABLE IV. Scaling function parameters 



Parameter / System Percolation (pj) Percolation (m) Random Partitions Au + C Scaled Au + C 

al 0.8±0.2 1.8±0.2 0.75±0.5 5.9±0.1 1.9±0.2 

61 -1.0±0.1 -2.6±0.3 -3.0±0.5 -7.5±0.3 -3.1±0.3 
ci 0.5±0.2 1.8±0.2 4.0±0.5 3.2±0.1 1.8±0.2 
a 2 1.0±0.1 0.3±0.03 0.75±0.5 0.8±0.2 6.9±0.7 

62 -0.5±0.1 O.liO.01 -3.0±0.5 -1.2±0.4 -ll.Oil.O 
c 2 0.8±0.2 1.1±0.1 4.0±0.5 1.5±0.2 3.8±0.4 
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FIG. 20. Results of the scaling function analysis with the 
critical point equal to half the value as determined from the 
Fisher r-power law analysis for: (a) percolation as a function 
of z(pi), (b) percolation as a function of z(m), (c) random 
partitions as a function of z(m) and (d) Au + C multifrag- 
mentation as a function of z(m). 



FIG. 21. Results of the scaling function analysis with the 
critical point equal to twice the value as determined from the 
Fisher r-power law analysis for: (a) percolation as a function 
of z(pi), (b) percolation as a function of z{m), (c) random 
partitions as a function of z(m) and (d) Au + C multifrag- 
mentation as a function of z(m). 



mined here is the scaling function for percolation in three 
dimensions, i.e. it is universal for three dimensional per- 
colation independent of size. The scaling functions for 
pi and m determined above agree well with the scaled 
cluster distributions of different size lattices, see Figure 
23, and can be used to predict the behavior of the second 
moment for any size lattice |fn| . In the same spirit, the 
scaling function determined here for gold multifragmen- 
tation is the scaling function for charged nuclear matter 
which describes the cluster distributions produced in the 
multifragmentation of any nucleus, not just the excited 
gold remnant discussed in this work. With the knowledge 
of the form of the scaling function various other quanti- 
ties can be determined as illustrated in section III and 
shown below. 

The cluster distributions for the random partitions is 
fit, by eye, with the same empirical parameterization as 
in eq. ( |3^ ) see Figure 19c. The random partitions cannot 
be described by eq. (|33|). The solid curve in Figure 19c 
will be used in the following section to demonstrate the 
failure of the scaling analysis, as is also seen here, when 
applied to a system where a continuous phase transition 
is absent. 

Finally, a consistency check in this analysis is the 
agreement between the location of the peak in the scal- 



ing functions and the values of z max determined in the a 
analysis, see Table III. 

4- 'y-power law from the scaling function 

The behavior of kt or M2 can be derived from the 
functional form of the scaling function and the critical 
parameters via eq. Jl~7|). Performing the integration in 
eq. ([l7|) using the functional form of the scaling function 
determined above yields a direct calculation of the critical 
amplitudes, T± via eq. (p^). The critical exponent 7 
is calculated from the values of r and a via a scaling 
relation in eq. (|i"g|). Combining these two, T± and 7, 
it is possible to calculate the 7-power law that describes 
the behavior of the second moment. This calculated 7- 
power law can then be compared to the behavior of M2 
as measured from the cluster distribution. Figure 24a, 
b and d shows the agreement between the measured M% 
data (largest cluster omitted in the liquid region) and the 
calculated 7-power law curves for percolation (pi and m) 
and gold multifragmentation, respectively and Tables II 
and III list the results. 

The values of 7 determined via the scaling relation 
in eq. (|T^) for percolation (p; and m) show approximate 
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FIG. 22. A measure of the scatter in the scaling function 
analysis as a function of the choice of the critical point for: (a) 
percolation (pi), (b) percolation (m), (c) random partitions 
and (d) Au + C multifragmentation. 



FIG. 23. The scaling function analysis for cluster distribu- 
tion from a percolation lattice of Ao — 64 as a function of (a) 
Pi and (b) m. Solid curves show the fitted scaling function 
from the Ao — 216 lattice. 



agreement with the accepted value of 1.8. The high value 
of a extracted above leads to a low value of 7 here. Fig- 
ures 24a and b also show the behavior of the second mo- 
ment of a 250, 047 site lattice. The power law predicted 
using the scaling function determined with a 216 site lat- 
tice shows rough agreement with the measured M2 of the 
larger lattice in both the amplitude (r±) and exponent 
(7). There is approximate agreement between the pre- 
dicted power law and the measured M2 of the smaller 
lattice over some region in e that is neither too near to, 
nor too far from the critical point, e = 0. It is this region 
that will be determined, independently, in the following 
section. 

For the percolation (pi and m) system, the disagree- 
ment between the measured Mi data and the calculated 
curves is due to two well known reasons: far from the crit- 
ical point, the assumptions of scaling are no longer valid 
and the analytic background overwhelms the singular be- 
havior. Near the critical point finite size effects dominate 
M2, limiting the sizes of the large clusters which make the 
most significant contribution. In contrast, the r-powcr 
law was observed at the critical point because it is deter- 
mined by smaller clusters which suffer the least from the 
finite size effects. 

Figure 24c shows the results when this analysis was 
applied to random partitions. The power law predicted 
from the scaling function analysis applied to the cluster 



distribution of the random partitions fails to reflect the 
behavior of the measured second moment. This is not 
surprising as the random partitions presented here are 
not the result of a system undergoing a continuous phase 
transition. The disagreement observed in Figure 24c then 
serves as an indication of how this particular analysis 
probes for the presence of a continuous phase transition. 
This figure shows the results of this analysis for a system 
with no phase transition, while Figures 24a and b show 
the results of this analysis on a system where such a phase 
transition is present. 

The results of this analysis when applied to nuclear 
multifragmentation are shown in Figure 24d. In this case, 
the comparison to the predicted 7-power law is neither 
as good as that for percolation nor as poor as that for the 
random partitions. It shall be shown in section VI that 
considerable improvement can be achieved if account is 
taken of the changing system size, Aq (m) , and finite size 
scaling effects. 

The approximate agreement between the predicted 7- 
power law and the measured Mi behavior is in keeping 
with the behavior expected for small systems undergoing 
a continuous phase transition, e.g. the percolation sys- 
tem. The multifragmentation results are clearly different 
that then results of a system without a continuous phase 
transition, e.g. random partitions. 
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FIG. 24. Results of the 7-power law determined from the 
scaling function for: (a) percolation (pi), (b) percolation (m), 
(c) random partitions and (d) Au + C multifragmentation. 
Solid lines show the 7-power law predicted from the scaling 
function, a and r. Filled circles show the second moment as 
a function of e. Open squares in (a) and (b) show the second 
moment behavior of a percolation lattice of size Ao = 250, 047. 



5. ^-matching 

In the previous works the procedure for determining 
critical exponent values and the location of the criti- 
cal point from the cluster distribution was based on a 
method of matching exponent values on both sides of 
the critical point ||, |30|. The idea was to find the re- 
gion on either side of the critical point where the power 
law behavior predicted by the scaling function holds. As 
is seen in Figure 24 there is some intermediate e region 
where the second moment data are described by a power 
law, a region where the M 2 behavior is dominated by the 
7-power law and all other effects are small in compari- 
son. In earlier percolation studies J60[ general guidelines 
based on the correlation length and size of the fluctua- 
tions were used to find the boundaries in e of the regions 
to be fit. In nuclear multifragmentation analyses || it 
was impossible to use such guidelines. Instead a method 
was developed that searched for regions best fit by power 
laws and determined the location of the critical point 
and exponent values simultaneously. The values of the 
critical exponents and the normalizations associated with 
power laws were obtained from the best fit power laws in 
those regions. As with the previous analyses presented 



in this paper, this method of exponent matching does 
not select a particular value of a critical exponent or the 
critical point. Instead the values found are the outcome 
of an unbiased procedure. 

The method is as follows. A choice of the critical point, 
p c or m c was made. From this choice plots such as those 
shown in Figure 24 were made. Then fitting boundaries 
in e were chosen. The fitting range was defined by e£ ar 



and 



For example, on the gas side of the critical 



point a fit of /n(M 2 ) versus Zn(|e|) was made for all data 
with \e\ ear \ < \e\ < \e f + ar \. The slope of the resulting 
linear fit was recorded as 7+, the offset as ln(T + ) and 
the goodness of fit as xl+- The same procedure was 
applied to the liquid side of the chosen critical point, 
recording 7_, in(r_) and xt_. For each choice of the 

critical point, several choices of fitting regions, e£ ar and 
£ near were m ade and results recorded. Five parameters 
were chosen for each set of power law regions examined: 
e£ ar , e^f ar and p c or m c . 

The 7-power law fit regions and critical point loca- 
tions were evaluated by demanding that: (1) they yield 
7+ and 7_ values that matched each other to within the 
error bars on those values returned by the fitting routine 
and (2) that the xt °f the fits were in the lowest quarter 
of the distribution resulting from all the fits which satisfy 
condition (1). The results from the power law fit regions 
that passed these two criteria were then histogrammed 
and average values for all quantities concerned were de- 
termined. The results are summarized in Tables I, II and 
III and shown in Figure 25. 

The lines plotted in Figure 25 do not result from any 
single fit, but display the average results for j± and r± 
that have satisfied conditions (1) and (2). The points in 
Figure 25 are the measured second moment for the par- 
ticular cluster distribution in questions plotted against e, 
which depends on the average value of p c or m c that sat- 
isfies conditions (1) and (2). Therefore the lines in Figure 
25 should not be interpreted as a fit to the data points 
shown in the same figure, but as the average results from 
the 7-matching procedure. Full circles in Figure 25 show 
the average fitting regions that satisfy conditions (1) and 
(2). 

For percolation pi the value of 7 determined in this 
manner is within a few percent of the value determined 
in and the infinite lattice value. The ratio of T + /r_ 
determined by this method, a ratio that depends on the 
universality class of the system in question, is also in 
agreement with the infinite lattice value and the T± val- 
ues predicted by the scaling function, see Table III. The 
value of p c determined here is within 15% of the value de- 
termined in a previous analysis of the L = 6 lattice ]6C| ] 
and the value determined above in the Fisher r-power 
law analysis, see Table II and Figure 25a. 

The results for the analysis of percolation with m as 
a measure of the control parameter are worse that the 
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FIG. 25. Results of the 7-matching analysis for: (a) per- 
colation (pi), (b) percolation (m), (c) random partitions and 
(d) Au + C multifragmentation. Solid lines show the average 
7-power fit returned by the procedure. Open circles show the 
entire second moment behavior, filled circles show the average 
fitted region returned by the procedure. See text for details. 



results when the natural control parameter pi is used, the 
difference in 7+ and 7_ was: Aj m = 0.06±0.1 compared 
to Aj Pl = 0.0 ± 0.3. This is to be expected because for 
each value of pi there is some spread in the resulting 
values of m, so that binning in m groups together events 
with different values of pi. There is also a non-linear 
relation between the average values of p\ and m |7l|| . In 
spite of these two effects the results of the analysis in 
section IV-B-4 suggests that vestiges of the signature of 
a phase transition are still present even when m is used 
as the control parameter. That is also the case in the 
present analysis. Table II shows that the 7 value agrees, 
within error bars, with the infinite lattice value. The 
values of the critical amplitudes, T±, do not yield a ratio 
that agrees with the infinite lattice value. This is due to 
the non-linear mapping of pi onto m and is discussed in 

When the 7-matching procedure was applied to the 
cluster distributions from random partitions a very lim- 
ited amount of trial fits passed the combined tests of (1) 
and (2). The results compared poorly to the percola- 
tion results. At best the values of 7+ and 7_ match to 
within 20% of the average value of 7, compared to perfect 
matching for percolation p and matching within 5% for 
percolation m. The value of the critical point, m c , re- 



turned from this analysis also compared poorly to other 
outcome of previous analyses, see Table I. Finally, while 
fit regions for all systems were limited, the fit regions are 
the smallest and the poorest of quality for the random 
partitions. 

The results of the 7-matching analysis applied to mul- 
tifragmentation data has been published in rcf. In 
that work the data were contaminated by the inclusion 
of prompt nucleons; prompt nucleons are excluded from 
consideration in this work. In that work the second mo- 
ment of the cluster distribution was determined based on 
the charge of a cluster rather than its mass as is done in 
this work. Previously, the second moment was generated 
from a cluster distribution that was not normalized to the 
changing size of the system as is done here. Furthermore 
the prior analysis consisted of only one quarter of the to- 
tal number of events used in the present analysis. Thus 
the current analysis has higher statistics, has been freed 
of prompt nucleons, has a second moment that has been 
constructed with the masses from the cluster distribution 
and a cluster distribution that has been normalized to the 
changing system size. The exclusion of prompt nucleons 
and normalization to the changing system size are an ef- 
fort to address the criticisms raised in and rebutted 
in 74 1. When the 7-matching procedure was applied to 



the data presented in this paper essentially the same re- 
sults as presented in ref. || were recovered. See Table II 
and Figure 25d. One difference observed is in the value 



of the critical point returned, 



,94 _ 



26±1 reported in 



ref. |8j] and = 21±2 reported in this work. The dif- 
ference is not as great as it appears to be. The origin 
of the published value of m 9 lies in picking the peak of 
the distribution of m c values that satisfied conditions (1) 
and (2) as the location of the critical point. The value 
was estimated based on the height of the peak and the 
error based on the width of the peak. The mean and 
RMS of the m c distribution in ref. || suggest a value of 
the critical point of m 9 c A — 25±3. This value agrees, to 
within error bars, with the value of m 99 presented here. 
The relatively small shift in m c can then be understood 
to arise from the differences in the data sets. Noting 
this it is clear that the present 7-matching analysis is in 
agreement with the previous work. 

The results of the present work are, again, in keeping 
with the expected results of a small system undergoing a 
continuous phase transition. There is some region where 
matching 7 values can be obtained, some regions in e 
where th 7-power law overwhelms all other effects. The 
fits in Figure 25d are of quality than those for random 
partitions in Figure 25c and cover a greater range. When 
compared to the percolation m results the multifragmen- 
tation data compare favorably in terms of overall good- 
ness of fits, width of fit region and matching of j± . See 
Table II. The location of the critical point returned by 
this analysis also compares well with the location from 
other analyses. See Table I. 
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V. CORRECTIONS TO SCALING 

In the last section it was seen that the 7-power law 
and the data for the second moment in all systems have 
agreed over only a limited area. To some degree this is to 
be expected. Near the critical point, assumptions valid 
for thermodynamic systems are invalid for the finite sys- 
tems discussed in this work. For that reason, finite size 
effects dominate at the critical point and the second mo- 
ment merely peaks instead of diverging. Far from the 
critical point other effects come into play. The scaling 
assumptions inherent in the FDM are valid only in the 
neighborhood of the critical point. The size of this neigh- 
borhood is somewhat ill defined and seems to depend on 
many factors, e.g. the quantity in question, the nature 
of the system, the size of the system and so on. Scaling 
behavior in physical systems can be observed over a wide 
range in temperatures and densities. This is most ele- 
gantly illustrated in the Guggenheim Plot jn| of scaled 
temperature (T/T c ) as function of scaled density {p/p c ) 
for several different gases (Ne, Ar, Kr, Xe, N2, O2, CO 
and CH4). In that plot the data collapse onto a curve 
that is well described by a power law with an exponent 
of (3 = 1/3. The range in validity of this agreement be- 
tween data and power law is shown on the Guggenheim 
Plot to be over a range of AT ~ 0.5T C and Ap ~ 2.5p c . 
However, another system, the combination of isobutyric 
acid and water, shows the Guggenheim type of scaling 
only very near the critical point |76j ] . Already when the 
range considered is AT ~ 0.04T C and Ap ~ 0.0fp c cor- 
rections to scaling can be observed. To that end, higher 
order corrections to scaling are now examined in order to 
determine if fits such as those shown in previous sections 
can be improved. However, any improvement comes at 
the expense of more fit parameters and assumptions. 

To fully explore corrections to scaling in the context of 
the present systems where the cluster distributions serve 
as the main observable the FDM is revisited in a fashion 
employed in references Q and |fT7[ . Assuming coexis- 
tence eq. (||) is then re- written as 



n Af (e) = q Aj r (/<>(*) + Aj n fi{z) + . . .) 



(34) 



where fi(z) is the correction-to-scaling function and f2 is 
the correction-to-scaling exponent. The form of eq. j3^ ) 
anticipates the presence of a second function of z. In sec- 
tion IV-3 it was found empirically that both the scaled 
percolation and multifragmentation cluster distributions 
(uA f (e)/qoAj T ) could be reasonably well described by 
the sum of two gaussians, eq. (|33|). In that treatment, 
the amplitude of each gaussian was a constant, a,i and 
ai. If Af is restricted to a single value, the prescrip- 
tion give by eq. ([54]) is equivalent to that of eq. (33). 
Eq. ([mJ) predicts that there should be an ordering to 
the scaled cluster distributions, i.e. smaller cluster sizes 
should lie above the larger clusters due to the correction 



term. This can be observed in Figures I9a and b in the 
neighborhood of the maximum of the scaled data for the 
percolation systems. In the tails of the distribution, ei- 
ther large or small cluster production is suppressed. In 
the case of multifragmentation data, Figure 19d, the or- 
dering is generally observed where the statistics are ade- 
quate, namely, near z = 0. The ordering of the random 
partitions implies that < 0. 

From eq. ([54]) it possible to derive the corrected 
isothermal compressibility (second moment) power law. 
Following the method in section III leads to: 
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(35) 

which is usually simplified and written as 

KT ~r ± | e p(l + a± | e | A ). (36) 

Now the overall amplitude, r± is given by the first inte- 
gral, and the correction-to-scaling amplitude is given by 
the second integral divided by the first. The correction- 
to-scaling exponent is A = Q/a. 

Using eq. ( |36| ) to fit the second moment distribution 
would lead to determining four fit parameters: two am- 
plitudes and two exponents. To explore the effects of 
corrections to scaling an assumption was made as to the 
universality class of the system in question and thus the 
choice of exponent values. For the three dimensional per- 
colation universality class A = 1.22 [ p8[ , J77J and for the 
three dimensional Ising universality class A = 0.56 |78||- 
pOfl . The amplitudes were left as free parameters and 
the second moment of the cluster distributions were fit. 
The value of the critical point determined from the 7- 
matching analysis of section IV-B-5 was used for each 
system. Figure 26 shows the results. 

For the percolation system with the corrections-to- 
scaling a better fit to the second moment data was pos- 
sible over a greater range in e, up to twice the range of 
the average fitted region in the 7-matching analysis, see 
Figures 26a and b. The fits still failed to reproduce the 
behavior of M<z near the critical point where finite size 
effects dominate the system. Table III lists the results 
for the critical amplitudes, T±. The agreement in the 
critical amplitudes determined in this analysis and the 
amplitudes from the 7-matching analysis is due to the 
agreement of the behavior of eq. ([36]) and the 7-power 
law from the 7-matching analysis over the region in e de- 
termined by 7-matching. Thus the 7-matching analysis 
finds regions that are the least affected by higher-order 
corrections to scaling. 
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FIG. 26. Results of the correction-to-scaling analysis for: 
(a) percolation (pi), (b) percolation (m), (c) random parti- 
tions and (d) Au + C multifragmentation. Open circles show 
the entire second moment behavior, filled circles show the fit- 
ted region used in this procedure. See text for details. Dotted 
lines show the 7-matching power law. Dashed (solid) lines 
show the resulting correction-to-scaling 7-power law for the 
3D percolation (3D Ising) universality class. 



For the random partitions an improvement is only ob- 
served for the high multiplicity region where a better 
fit over a larger range was obtained for both choices of 
universality class, see Figure 26c. The low multiplicity 
events showed no such improvement partly due to the 
limited range in e available. Both the three dimensional 
percolation and three dimensional Ising exponents were 
used in this analysis for the random partitions and the 
multifragmentation data. Both choices of universality 
classes showed similar results. The lack of effect of cor- 
rections to scaling is to be expected in a system that does 
not follow FDM like scaling laws. 

The multifragmentation data also showed improve- 
ment resulting in a better agreement between the fits 
and the M2 data points over a larger range in e, see Fig- 
ure 26d. The improvement was observed for both choices 
of universality classes thus indicating this analysis is in- 
sensitive to the differences 
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82 1, though the good- 



ness of fit was better for the choice of the three dimen- 
sional Ising exponents over the same fit regions (3D Ising: 
0.4 and xl- = 0.7; 3D percolation: xl+ = 0-7 and 
1.1). At this level of analysis it appears that cor- 
rections to scaling improves the fits for the 7-power law. 
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Whether this is due to the presence of a continuous phase 
transition in nuclear multifragmentation, or merely extra 
terms in the fitting function remains an open question. 



VI. CHANGING SYSTEM SIZE AND FINITE 
SIZE CORRECTIONS 

It has been pointed out that the previous statistical 
analysis of gold multifragmentation B ignored the chang- 
ing size of the system ]73| ] . To first order this may have 
been a reasonable procedure Jt3] as many statistical sig- 
natures of a continuous phase transition have been ob- 
served both before and after the scaling to account for 
the changing system size has been performed; e.g. the 7 
power law shown here and in previous works agree well. 
However, the data collapse in Figure 19d is qualitatively 
not as great as that shown by the percolation (m) system 
in Figure 19b and the agreement between the calculated 
and measured M2 behavior in gold multifragmentation 
in Figure 24d is qualitatively not as good as that shown 
by percolation (m) in Figure 24b. In this section, the 
effect of the changing size of the system is explored and 
accounted for. 

The size of the multifragmenting system is shown in 
Figure 27a after ref. pj] |. An approximately linear re- 
lation between the system size, Aq, and m was found: 
Aq = 199 — 1.6m, see Figure 27a. The functional form 
of Ao(m) was used in the following analysis to account 
for the changing system size in an average way, i.e. not 
on an event-to-event basis. 

If the multifragmenting system is assumed to be a sys- 
tem undergoing a phase transition, then the theory of 
finite size scaling of the critical point |33|- |^6| asserts 
that the effective critical temperature, T c (Aq), changes 
as a function of the system size. Coupling this with the 
changing size of the system indicates that at each value 
of to the value of T^^(Ao) is different. 

The value of T^'* (Aq) can be determined in the follow- 
ing manner. First a relation between the multifragment- 
ing system's temperature, T , and m must be determined. 
Again from ref. |5l[| a relation can be found, see Figure 
27b. Two estimates of T were made, one from a Fermi 
gas (uncorrected for the effect of expansion energy), T, 
and the other from an isotopic yield ratio thermometer, 
Tf, |87j| . These temperatures give an approximate indi- 
cation of the initial and final temperatures of the sys- 
tem. Figure 27b shows that Tj is well fit by a quadratic 
function, while Tf is well fit by a linear function. An- 
other linear function reproduces the average of T and Tf. 
T = 3.0 + 0.14to, this was used for the following analysis. 
The critical temperature of infinite nuclear matter was 
assumed to be T c °° = 17 ± 1 MeV §§]. 

From the Fisher r-power law analysis the value of the 
multiplicity at the critical point was determined to be 
to = m c — 22 ± 1. The system size at that point is then 
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FIG. 27. (a) Behavior of the size of the system, Ao, with 
respect to the total charged particle multiplicity, m. The solid 
line is a fit to the changing system size, (b) measures of the 
temperature, T, as a function of m: □ Ti from a Fermi gas and 
O Tf from an isotope thermometer. The solid curve shows a 
fit to Ti, dotted curve a fit to Tf and the dashed curve a fit 
to(T> = (T<+2»/2. 



A (m = 22) = 164 ± 2 and the temperature is T = 6 ± 2 
MeV. This indicates that the critical temperature for a 
nuclear system with 164 nucleons, T^-f(A ), is approxi- 
mately 6 ± 2 MeV. 

According to theory, to first order the critical point 
scales with system size as: 



(T c e ^(A )-T c °°)/T c c 
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(37) 



where d is the Euclidian dimension of the system and v is 
the so-called hyperscaling exponent. At the smallest of 
system sizes higher order correction terms may play an 
important role in the scaling of the critical point J|?| . It 
is assumed in this procedure that the above mentioned 
finite size scaling dominates all other effects, including 
those which arise due to the charge of the protons present 
in the excited remnant. This is in keeping with the gen- 
eral philosophy of the theory of critical phenomena where 
near the critical point the precise details of the system 
are irrelevant. 

If it is assumed that gold multifragmentation is the 
result of a continuous phase transition and that transi- 
tion occurs in three dimensions, d = 3, then using the 
hyperscaling relation |32| 



with the extracted values of a and r, gives v ~ 0.63±0.07. 
The coefficient b in eq. ( |37| ) can be determined using 
T^f{A Q {m c )), T c °°, d and v; resulting in b = -9 ± 2. 
Note that this value of b suggests that T*'* (Ao) = for 
Ao = 70 ± 40. This is a result of the form of eq. ( |3^ ) and 
the notion that the critical temperature lowers as the 
size of the system decreases. Presumably higher order 
effects not taken into account in eq. (37) will affect the 



location in system size, Aq, where the effective critical 
temperature vanishes. For the form of finite size scaling 
corrections shown in eq. (|37|), only b = — 1 yields an 
effective critical temperature that vanishes at Aq = 1. 

Now eq. ( |37| ) can be used to solve for T c e ^ (Ao(m)) 
and thus determine the scaled control parameter: 

e scaled = (Tff(A (m)) - T)/Tff(A (m)). (39) 

The analysis to extract the exponent a was performed 
with this e scaled by finding the peak in Af sized cluster 
production as a function of nA f (e scaled ) versus e scaled . 
Previously, it was argued that the largest cluster should 
be excluded for all values of e to account for finite size 
effects in the analysis to extract a. The analysis in this 
section directly accounts for finite size effects, thus the 
standard FDM formalism with respect to the largest clus- 
ter is followed. This results in a value of a — 0.65 ± 0.07 
and Zma l x d = — 11 ± 2, Figure 28a shows the resulting 
power law. 

The scaling function for gold multifragmentation was 
then plotted using the above corrections for the changing 
system size and finite size scaling, see Figure 28b. The 
data collapse is qualitatively better than in Figure 19d. 
The two gaussian parameterization of j( z scaled ) was fit 
to the scaled scaling function and is shown in Figure 28b 
with parameters listed in Table IV. 

Using the fitted parameterization of the scaling func- 
tion and other quantities, the 7-power law can be deter- 
mined as before. Figures 28b and c show the measured 
M2 of gold multifragmentation plotted as a function of 
A 7-power law was plotted on, not fitted, Figures 
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28b and c with 7 = 1.3±0.2 (from r, a and eq. (18)) and 
offsets determined via the scaled scaling function. This 
7-power law agrees with the measured M2 over nearly all 
of the e sca ' ed -range, the exception being near e scaled = Q 
where finite size of the system limits the maximum of 
M 2 . 

As stated above, this analysis makes no attempt to 
account for the Coulomb energy of the system in the for- 
mulation of finite size scaling. The effect of the Coulomb 
energy may evidence itself in the degree of collapse of the 
data. The fact that, to a large extent, the data do col- 
lapse without any explicit adjustment to the theory (i.e. 



formalation of the scaling variable, 



..scaled\ 



may then be 



an indication that the Coulomb energy is not a major 
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FIG. 28. Results of the system size / finite size scaling anal- 
ysis for Au + C multifragmentation. (a) The scaled a power 
law. (b) The scaling function, normalized cluster distribution 
versus z scaled \ The solid curve show the fitted scaling func- 
tion, (c) and (d) show the resulting 7-power law as predicted 
from the scaling function, a and r. Solid circles show the sec- 
ond moment plotted against the properly scaled e scaled an d 
solid lines show the predicted 7-power law. 



perturbation. Hence, it is perhaps justifiable to adhere 
to the Coulomb-free theory when scaling T c . 



VII. DISCUSSION AND CONCLUSIONS 

The focus in the present paper has been on cluster dis- 
tributions and the types of analyses which can shed light 
on their creation mechanism. In particular, attempts 
were made to identify procedures that can distinguish 
those distributions which are related to critical behavior 
from those which are not. While this question is easily 
answered for systems containing large numbers of con- 
stituents, it is much more difficult to address the case 
of interest here, namely systems with at most only a few 
hundred particles. In such systems, finite size effects play 
a large role. Therefore, two different computational mod- 
els, bond building percolation and a random partitions, 
have been used. It is well known that in the macroscopic 
limit, the former system possesses a continuous phase 
transition characterized by a unique scaling function and 
set of critical exponents while the latter system does not. 
In addition, data arising from the multifragmentation of 



gold nuclei has been studied using the same procedures. 
For this system, it is not known, a priori, whether a crit- 
ical point is present. 

Many cluster properties have been proposed as being 
suitable measures of critical behavior. Among these are 
the fluctuations in the size of largest fragment (Figure 1), 
peaking in the quantity 72 (Figure 3), peaking behavior 
in M 2 (Figure 5), Campi plots (Figure 6) and simple 
power law behavior in the cluster mass distribution for 
a particular value of the appropriate control parameter. 
It was seen that none of these measures, taken alone or 
together, was sufficient to distinguish a system possessing 
critical behavior from one which does not. 

The first procedure which produced different results for 
critical and noncritical systems was the single parameter 
power law fit to the cluster mass distribution (Figures 
12 through 15). For the percolation systems and for the 
multifragmentation data, it was shown that the one pa- 
rameter power law fit describes the data well only over a 
very limited range of the control parameter. The value 
of the control parameter where the power law fit is best 
is very close to where the abovementioned peaks occur. 
However, for the random partitions, this was not the case. 

If a system possesses a critical point, it is also expected 
to possess a scaling function that describes its behav- 
ior away from the critical point. Phase transition the- 
ory specifies how the argument of this function depends 
on cluster mass and distance from the critical point. If 
such a function exists, the theory permits the determi- 
nation of the critical exponent a. This determination 
was done for the percolation system yielding satisfactory 
agreement with its known value. The same procedure 
was applied to the random breakup model and to the 
multifragmentation data. It was clear from this anal- 
ysis that percolation and multifragmentation were very 
similar in many features, while the random system was 
distinctly different. See Figures 17 and 18. 

Again, if a system possesses a critical point, it is ex- 
pected to possess a scaling function that describes its 
behavior away from the critical point. Therefore, when 
the data is properly scaled, it should collapse onto a sin- 
gle curve. Figure 19 shows the amount of data collapse 
for the systems discussed here. The quality of the data 
collapse (Figures 19-21) reinforces the notion that the 
random breakup system is different from the others. Al- 
though the precise form of the scaling function is not 
dictated by phase transition theory, both the percola- 
tion system and multifragmentation data were satisfac- 
torily described by a sum of two gaussians. The random 
breakup model was not. 

The issue of finite size scaling was discussed. Unlike 
the other systems examined here, which had a fixed num- 
ber of constituents, the nuclear multifragmentation data 
originated from systems whose size varied monotonically 
with observed charged particle multiplicity. See Figure 
27a. Phase transition theory makes a prediction, eq. 
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( P7| ) , as to how the value of the control parameter at the 
critical point changes as a function of system size. Ap- 
plying eq. (H?]) produced an improvement in the quality 
of the data collapse for the scaling function and yielded 
a better prediction for the behavior of M^. See Figure 
28. 

Critical exponent values have been determined in an 
unbiased manner for each system. For both sets of anal- 
yses on the percolation clusters, the standard percolation 
exponents were recovered to within error bars. For the 
random partitioning, exponents could be extracted, but 
none that fulfilled well known scaling laws. The exponent 
values determined from the gold multifragmentation clus- 
ter distributions fulfill the scaling laws, to within error 
bars, and fall near the three dimensional Ising universal- 
ity class. 

The effect of secondary decays from hot initial frag- 
ments on the critical exponents has not been explicitly 
considered in this paper. In the SMM [p2| - fi5f such ef- 
fects become significant above E*/A = 7 MeV/nucleon. 
Thus r and 7, which are determined at lower excitation 
energies, will be unaffected in the SMM's fragment dis- 
tributions. The exponent a is determined by the multi- 
plicities at which individual light fragment yields attain 
their peak values. As shown in Figure 17d, the light- 
est fragments peak at large multiplicities, corresponding 
to excitation energies for which secondary decay are im- 
portant in the SMM. An SMM calculation Q indicates 
that the value of a from the SMM's fragment distribution 
was increased by about 70% due to this effect. However, 
it is unclear from that calculation that the effects of sec- 
ondary decay are as great in the experimental data as 
they are in the SMM. The SMM calculation over pre- 
dicts the yield of light fragments which could indicate 
that the SMM estimates of secondary decay are too se- 
vere. Corrections to the model independent quantities 
determined in this paper based on the SMM calculations 
are premature. See Appendix C for discussions of an- 
other set of model based interpretations. 

Although the multifragmentation data possess many 
of the gross and detailed characteristics that the percola- 
tion system does, it is not at all obvious why this should 
be the case. After all, real nuclei obey quantum mechan- 
ics, have varying binding energies per particle, and, most 
significantly, are charged. On the other hand, it is well 
known that near a critical point the details of the inter- 
action become unimportant and only the dimensionality 
of the system and the dimension of the order parameter 
are important. As noted in the introduction, the attrac- 
tive nuclear force bears a similarity to a van der Waals 
force. However, the Coulomb force is a long range force 
and imposes a natural limit to the size of stable nuclei. 
Thus, it is not clear to what extent a finite charged sys- 
tem can exhibit critical behavior when the macroscopic 
system cannot exist. The exact role of the Coulomb force 
in physical systems undergoing a change of phase is cur- 



rently of great interest Q, M and is, at this point, an 
open question. The philosophy of this paper has been 
to make use of phase transition theory as it applies to 
uncharged systems. What results for the analysis of the 
gold multifragmentation data bears great similarity to 
the results of the same procedures applied to a system 
known to possess a critical point. It is tempting then 
to conclude that multifragmentation is related to critical 
behavior occurring in a finite nuclear system. 



APPENDIX A: DISCUSSION OF POWER LAW 
IN RANDOM PARTITIONS 

As a demonstration of the power of the sort of scal- 
ing analysis presented above it is shown that the random 
partitions follow a simple power law of Na } ~ Aj . Fig- 
ure Al shows the scaled cluster distribution as a func- 
tion of m from clusters with Af > 3. The data nearly 
collapses to unity along the horizontal axis over the mul- 
tiplicity range for m > 5. The deviations are due to the 
constraints of m and finite size. With a simple scaling 
analysis the underlying power law describing the cluster 
distribution becomes clear. 
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FIG. 29. Results of a scaling analysis performed on the 
cluster distribution of random partitions. The scaled cluster 
distribution (nearly) collapses to unity over the entire multi- 
plicity range. 
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APPENDIX B: RIEMANN £ FUNCTION 
SUMMATION 



The value of qo used in this work based on the £- 
function was generated with VAX FORTRAN code using 
double precision and letting the sum run from 1 to 10 6 . 
The sum was terminated at this point in order to keep 
computing times within reason. For a value of r = 2.18 
summing to 10 6 gives a value of qo that is within 10% of 
the value when the sum is terminated at 10 10 , see Figure 
A2. Increasing the upper limit of the summation in the 
(■-function causes no significant changes in the analysis 
presented in this work. 



APPENDIX C: LATTICE GAS MODEL 
INTERPRETATIONS 

In a recent paper |Q, Gulminelli and Chomaz (GC) 
report on calculations made using a lattice gas model. 
Their results can be divided into two different categories: 
(1) an analysis of the thermodynamic quantities of the 
system; (2) an analysis of the statistical properties of the 
droplet or cluster distribution. 

In their paper GC performed canonical calculations 
on a three dimensional lattice of n sites (n = 216 or 
512) with periodic boundary conditions at a given tem- 
perature, T, and density, p. The density was varied by 
changing the number of particles in the system, Ao, with 
p = Aq/ti. For a choice of T, p was varied and the sys- 
tem's chemical potential, pi, was determined. At some 
values of T the p vs. p (p vs. ^4o) isotherm exhibited a 
back bend. As T increased the back bend disappeared. 
This was taken as evidence that a first order phase tran- 
sition (back bending p vs. ^4o) had culminated in a con- 
tinuous phase transition (flat p vs. Ao). The critical 
point, (T c , p c ), was thus defined. Via a Maxwell con- 
struction, GC also determined the boundary of the coex- 
istence zone. 

The theoretical ideas behind the procedure followed 
by GC are sound (with the exception of the questionable 
definition of volume in a system with periodic boundary 
conditions), however the execution is flawed. The error 
in execution lies in the method by which p was varied. 
GC varied p by holding n, the volume of the system, fixed 
and varying Aq, the number of constituents or "size" of 
the system. Therefore, at each value of p GC have a sys- 
tem of different size. It is well known that the effective T c 
of a system varies with its system's size fflj . Thus 
each time GC varied Aq to vary p they were dealing with 
a system with a different T c . With this understanding 
it is clear that a plot of p versus Ao as shown by GC in 
their Figure, la, which serves as the basis for their de- 
termination of the critical point, is difficult to interpret. 
T c cannot be deduced in the manner followed by GC. 
Similarly, it is impossible to determine p c for a system of 



Ao = 100 using information from a system of Aq = 200 
without accounting for the finite size scaling of the criti- 
cal point. This makes the thermodynamically extracted 
values of the critical point suspect and calls into question 
the validity of the coexistence line shown by GC in the 
lower part of their Figure 3. 

Via an analysis of the droplet distributions GC deter- 
mined the value of the critical exponent t — T max based 
on the peak in the production of different size droplets. 
The value of T c was determined coupling the T max results 
with a two parameter power law fit the droplet distribu- 
tion at various temperatures, r(T), at T c T max = t(T c ). 
The use of two parameter fits is generally an improper 
method to determine r 71 1- [Bsl. In this analysis GC 
determine T c on a system-by-system basis, avoiding the 
error of grouping systems of different size together. GC 
also extracted the exponent a and then used the above 
results to collapse the droplet data onto a single curve 
illustrating the scaling function that modifies the power 
law away from the critical point. While the data col- 
lapses onto a single curve, the value of the curve at T c 
is problematic. From GC's work the value of the scal- 
ing function at T c is f(T c ) ~ 0.5 which contradicts the 
r values extracted by GC according to the Riemann (- 
function relation between r and f(T c ) ]5S| ]; e.g. r ~ 2.2 
gives /(T c ) - 0.2, while GC's f(T c ) ~ 0.5 gives r - 2.7. 

Finally GC combine the results of their analyses and 
assert that the Kertesz line of r power laws continues into 
the coexistence region of the lattice gas and that this be- 
havior is observed due to the finite size of the system and 
is not observed in systems of much greater size. How- 
ever, due to the problems in GC's method of varying p 
to find the coexistence curve, such a claim is premature. 
Instead, their estimates of T c from the droplet analysis, 
which are shown as a string of points in their Figure 3, 
may reflect the finite size scaling of the critical point and 
not the Kertesz line extending into the coexistence zone. 
That each different size system appears to be at p c as well 
as T c , illustrated by the data collapse, may be explained 
by the relation p — p coex ~ (p — p c ) S ■ For a 3D Ising 
system 6 ~ 4.8. Thus, systems of finite size near their 
critical density are nearly on the coexistence curve and 
show the sort of data collapse illustrated in GC's Figure 
2§ ; 
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FIG. 30. Results of the summation to determine the offset 
go as a function of the upper limit of the summation. 
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